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Chapter  1 


Introduction 


1.1  Unsteady  Aerodynamics 


The  study  of  truly  unsteady,  high-excursion  and  high-Reynolds  Number  separated  flows 
over  submarines,  aircraft  dr  missiles  has  become  of  great  importance  in  the  analysis 
and  the  improvement  of  the  dynamic  performance.  Because  of  highly  complex,  three- 
dimensional,  turbulent  and  separated  nature,  standard  stability  derivative  techniques  fail 
to  capture  the  nonlinearities  in  these  flows  and  Computational  Fluid  Dynamics  (CFD) 
techniques  need  physical  models  that  can  resolve  the  complexities  of  such  flow  fields  in 
order  to  get  accurate  and  more  reliable  results.  Suitable  simulation  of  the  time-dependent 
maneuvers  in  the  wind  tunnels  is  not  only  important  for  understanding  the  physics  of 
complex  flow  phenomena,  but  also  supplies  the  necessary  information  required  for  devel¬ 
oping  the  realistic  unsteady  physical  flow  models.  The  DyPPiR  (Dynamic  plunge-pitch- 
roll)  model  mount,  a  computer  controlled,  three  degrees  of  freedom  robotic  arm  at  the 
Stability  Wind  Tunnel  of  Virginia  Tech,  provides  the  unique  capability  of  performing 
pre-programmed  general,  high-excursion,  large  scale,  high-Reynolds  Number  unsteady 
maneuvers  [4]. 

Dynamic  testing  has  been  an  important  part  of  design  and  validation  of  various  types  of 
craft  for  decades.  Typically  these  techniques  are  only  quasi-steady,  relying  on  very  small 
amplitude  sinusoidal  oscillations  that  can  describe  small-excursion  maneuvers  reasonably 
well  [2].  As  discussed  by  Wetzel  and  Simpson  [3],  there  is  a  significant  difference  between 
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quasi-steady  and  unsteady  aerodynamics.  In  a  quasi-steady  approach,  the  aerodynamics 
of  a  maneuvering  body  are  dependent  only  on  the  instantaneous  state  of  the  model  (cc 
angle  of  attack,  0  sideslip  angle,  control  surface  deflections,  etc.),  whereas  in  fully  general 
unsteady  aerodynamics,  explicit  time  dependency,  or  history  effects  are  also  included. 
Mathematically  the  distinction  between  quasi-steady  and  unsteady  aerodynamics  can  be 
shown  as  follows:  steady,  F(a,/3, ...);  quasi-steady,  F[a(t^),/3(t ), ...];  and  fully  unsteady 

where  F  and  G  can  be  a  dominant  flow  feature  such  as 
separation  location.  Here  f  is  the  non-dimensional  time  defined  by  Etkin  [5]: 


tref 


tUpo 

L 


(1.1) 


tref  represents  the  time  for  the  flow  to  pass  over  a  model:  L/Uoo-  The  non-dimensional 
time  t’  relates  the  unsteady  wind  tunnel  tests  to  the  real-time  maneuvers.  This  non- 
dimensionalization  follows  the  approach  that  is  used  to  non-dimensionalize  the  variables 
in  the  equations  of  motion  of  aircraft  dynamics  [5]  and  it  does  not  include  any  viscous 
terms.  However  the  onset  and  the  propagation  of  the  three-dimensional  unsteady  flow 
separation,  which  is  a  frequent  flow  phenomena  seen  over  a  vehicle  undergoing  a  typi¬ 
cal,  high-excursion,  and  truly  unsteady  maneuver,  are  strongly  influenced  by  the  viscous 
effects.  Reynolds  number  based  on  the  model  length  (or  any  other  appropriate  charac¬ 
teristic  length)  Rei  =  pUooL/iJ,  representing  the  ratio  of  the  inertial  forces  to  the  viscous 
forces,  is  another  important  parameter  to  be  considered  during  the  unsteady  wind  tunnel 
testing.  The  type  of  the  flow  separation  (laminar  or  turbulent),  and  the  transition  loca¬ 
tion  on  the  model  surface  are  strong  functions  of  the  Reynolds  number  until  a  critical 
value  is  reached.  Beyond  this  critical  value  of  Rcl,  the  separation  type  and  locations 
become  less  sensitive  to  the  Reynolds  number  effects.  In  some  experiments,  it  may  be 
impossible  to  simulate  the  real  flow  Reynolds  number  (such  as  Rcl  for  submarines,  big 
transport  airplanes  etc.)  due  the  model  size  or  the  wind  tunnel  speed  limitations.  In  this 
study  and  some  of  the  previous  unsteady  experiments  done  at  Virginia  Tech,  specially 
designed  boundary  layer  trips  are  used  to  simulate  higher  Reynolds  number  flows  and  to 
fix  the  transition  location  that  will  produce  a  less  Reynolds-number-sensitive  separation. 

Besides  quantitative  differences  between  the  steady  and  the  unsteady  aerodynamic  pa¬ 
rameters  such  as  the  skin-friction,  force  and  moment  coefficients,  the  qualitative  nature 
of  the  unsteady  flow  fields  also  differ  from  the  steady  ones.  In  fact,  the  quantitative  varia¬ 
tions  are  the  results  that  originate  from  flow  topology  differences  between  the  steady  and 
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the  unsteady  flow  flelds.  As  an  example,  in  the  crossflow  separation  phenomena  discussed 
in  this  study  and  previous  prolate  spheroid  work  [2],  the  history  effects  mainly  cause  time 
lags  between  the  unsteady  and  the  corresponding  quasi-steady  flow  fields,  which  can  be 
quantified  by  measuring  the  flow  separation  locations.  Approximating  these  time  lags 
by  the  models  that  can  explain  the  real  physics  of  the  three-dimensional,  turbulent,  un¬ 
steady  flow  separation  is  a  big  challenge  in  the  field  of  aerodynamics  as  well  as  acquiring 
unsteady  data.  One  of  the  main  purposes  of  this  study  is  to  apply  such  an  approximation 
to  the  unsteady  separation  locations  measured  on  the  DARPA2  model. 


1.2  Three-Dimensional  Separations 

Both  the  steady  and  the  unsteady  flow  fields  over  the  DARPA2  model  used  in  this  study 
are  dominated  by  the  three-dimensional  flow  separation  after  a  certain  angle  of  attack. 
As  described  by  Simpson  [6],  separation  is  the  entire  process  of  departure  or  breakaway, 
or  the  breakdown  of  the  boundary  layer  flow.  An  abrupt  thickening  of  the  rotational-flow 
region  next  to  a  wall  and  significant  values  of  the  normal-to-wall  velocity  component  V 
must  accompany  breakaway,  or  otherwise  this  region  will  not  have  significant  interaction 
with  the  inviscid  free-stream  flow.  As  can  be  seen  from  figure  1.1,  surface  skin-friction 
lines  converge  on  each  side  of  the  separation  line.  Along  this  separation  line,  there  exists 
a  stream  surface  across  which  no  flow  from  one  side  of  the  separation  line  can  pass  to 
the  other  side. 

Three-dimensional  (3-D)  separations  can  be  classified  in  three  groups  with  respect  to  their 
topology  and  kinematics  (figure  1.1):  (a)  horseshoe  type;  (b)  Werle  type;  and  (c)  crossflow 
separations.  In  Yates  and  Chapman  [1],  horseshoe  and  Werle  type  separations  are  defined 
as  global  separations  (or  closed  separations).  These  always  have  a  saddle  point  of 
separation  on  the  surface.  Figure  1.1  (c)  shows  the  local  or  the  crossflow  separation 
topology.  No  3-D  critical  points  on  the  surface  or  in  the  flow  can  be  observed  in  such 
type  of  flows  and  the  center  of  a  separated  vortex  structure  has  minimum  streamline 

curvature. 

The  horseshoe  type  and  the  crossflow  separation  are  of  particular  interest  in  this  study: 
the  flow  field  in  the  vicinity  and  the  downstream  of  the  sail  is  dominated  by  the  horseshoe 
type  separation  and  the  crossflow  separation  is  the  main  flow  characteristic  on  the  leeside 
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of  the  DARPA2  model  at  an  angle  of  attack  for  steady  and  the  unsteady  measurements. 
The  circumferential  pressure  gradient  is  the  dominant  factor  in  the  crossflow  separation 
process.  Figure  1.2  shows  the  on-axis  cross  sectional  view  of  the  secondary  flow  stream¬ 
lines  from  a  crossflow  separation:  SI  is  the  primary  separation  location  and  S2  stands 
for  the  secondary  separation  location.  R1  and  R2  locate  the  primary  and  the  secondary 
reattachment  positions.  At  low  angles  of  attack,  only  the  primary  separation  is  observed 
on  the  leeward  side  of  the  model  and  as  the  angle  of  attack  increases  the  primary  separa¬ 
tion  line  moves  towards  the  windward  side.  After  a  certain  angle  of  attack,  the  secondary 
separation  line  on  the  leeward  side  of  the  model  can  be  observed. 

Besides  qualitative  verification  of  the  crossflow  separation,  one  is  also  interested  in  finding 
the  actual  values  of  the  separation  locations.  Simpson  [6]  shows  that  the  wall  shear 
provides  lower  order  information  about  the  separation  location  than  the  pressure  by 
analyzing  the  continuity  and  the  momentum  equations  near  the  wall; 

Time-averaged  wall  pressure  measurements  are  relatively  insensitive  to  the  flow  separa¬ 
tion  and  only  massive  separations  are  detected  from  time-averaged  wall  pressure  mea¬ 
surements  [6].  In  order  to  use  equation  1.2  rigorously,  the  entire  wall  shear  direction 
and  the  magnitude  fields  should  be  known  to  perform  the  divergence  operation  on  the 
right  hand  side  for  finding  normal-to-wall  velocity  V  at  all  points  on  the  body  and  to 
determine  the  separation  location.  This  approach  is  impractical  due  to  the  large  number 
of  sensors  that  should  be  used.  As  an  alternative,  Simpson  et  al.  [7]  showed  that  the 
local  minimum  of  the  circumferential  skin-friction  distribution  is  a  good  approximation  of 
separation  location.  By  using  this  fact,  directionally  insensitive  hot-film  sensors  are  used 
to  measure  the  magnitudes  of  the  skin-friction  on  the  model  surface  in  this  study.  The 
local  minima  are  used  for  the  determination  of  the  steady  and  the  unsteady  separation 
locations. 

1.3  Previous  Studies 

Table  1.1  gives  a  summary  of  the  test  conditions  from  some  previous  unsteady  three- 
dimensional  aero/hydrodynamic  experiments  relevant  to  this  study.  In  the  last  row  of 
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Table  1.1:  Summary  of  the  test  conditions  from  some  unsteady  aero/hyrodynamic  experiments 
(modified  from  Wetzel  and  Simpson  [3]) 


Authors 

Model 

Tunnel 

Rei 

OL  —  OLyt 

b'oo 

a  range 

Measurements 

Gad-el-hak 

ogive 

water 

40,000 

1.05 

0  °  to  30  ° 

flow  vis. 

and  Ho  [8] 

cylinder 

Montividas 

cone 

wind 

56,650 

0.7 

0  °  to  90  ° 

flow  vis.. 

et  al.  [9] 

cylinder 

wake  LDV 

Smith  and 

Ogive 

wind 

1,200,000 

0.0405 

-15°  to  105° 

flow  vis.,  force 

Nunn  [10] 

cylinder 

&  moments 

Panzer  et 

hemisphere 

wind 

229,000 

0.0065 

15  °  to  30  ° 

wake  LDV 

al.  [11] 

cylinder 

Panzer  et 

hemisphere 

water 

75,000 

0.1 

15  °  to  30  ° 

wake  LDV 

al.  [11] 

cylinder 

Brandon  and 

F-18 

wind 

1,600,000 

0.0364 

0 

O 

00 

o 

o 

o 

t-H 

1 

flow  vis.,  force. 

Shah  [12] 

&  moments 

Wetzel  and 

6:1  prolate 

wind 

4,200,000 

0.047 

0  °  to  30  ° 

hot-film 

Simpson  [3] 

spheroid 

Whitfield  [13] 

DARPA2 

wind 

5,500,000 

0.076 

0  °  to  25  ° 

force  & 

Suboff 

moments 

Hosder 

DARPA2 

wind 

5,500,000 

0.071 

1°  to  27° 

hot-film 

(current  study) 

Suboff 

this  table,  information  about  the  current  work  is  presented.  The  maneuvers  performed 
for  all  the  experiments  shown  are  the  pitchup  motions.  The  model  geometry,  experi¬ 
mental  facility,  Reynolds  number  Rei,  non-dimensional  pitch  rate  a,  a  range,  and  the 
measurement  technique  of  each  study  can  be  compared  by  using  table  1.1.  In  terms  of 
Rei,  the  prolate  spheroid  study  of  Wetzel  and  Simpson  [3],  DARPA2  Sub  off  study  of 
Whitfield  [13]  and  the  present  study  can  be  considered  as  the  only  experiments  that  have 
high  Reynolds  number  well  above  the  critical  value.  Therefore  the  separation  locations 
in  these  experiments  are  least  likely  to  be  sensitive  to  the  Reynolds  number  effect.  These 
three  studies  are  performed  at  Virginia  Tech  Stability  Wind  Tunnel  with  the  DyPPiR 
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and  the  boundary  layer  tripping  are  used  for  all  three.  The  current  study  uses  the  same 
model  geometry  and  the  maneuver  type  as  the  ones  in  Whitfield’s  force  and  moment 
study  [13].  In  fact  these  two  experiments  are  complementary  in  the  sense  that  different 
measurement  techniques  are  used  to  describe  the  whole  picture  of  the  same  unsteady 
phenomena  over  the  same  geometry. 

Among  the  previous  studies,  the  unsteady  crossflow  separation  location  measurements  on 
a  maneuvering  6:1  prolate  spheroid  model  by  Wetzel  and  Simpson  [3]  can  be  considered  as 
the  most  similar  work  to  the  one  presented  here  in  terms  of  the  experimental  measurement 
technique  and  the  unsteady  high-excursion  maneuvers  performed.  In  that  study,  they 
have  determined  significant  lags  in  the  unsteady  flow  separation  locations  on  the  prolate 
spheroid  undergoing  pitch-up  and  turning  maneuvers  compared  with  the  steady  data. 
Wetzel  and  Simpson  [3]  also  found  that  a  first-order  lag  model  fits  the  unsteady  data. 
Present  work  and  the  prolate  spheroid  study  are  the  only  experiments  that  utilize  fine 
spatial  resolution  surface  hot-film  measurements  to  determine  the  unsteady  skin-friction 
magnitudes  and  the  separation  locations  on  a  maneuvering  body.  Most  of  the  previous 
work  uses  different  techniques  and  focuses  on  flow  features  other  than  separation,  such 
as  vortex  location  and  breakdown  [3] . 

1.4  Present  Work 

In  the  present  work,  unsteady  turbulent  surface  flow  on  a  maneuvering  DARPA2  subma¬ 
rine  model  is  studied.  Hot-film  sensors  are  used  to  measure  the  steady  and  the  unsteady 
sl^in-friction  magnitudes  over  the  body  surface.  Local  minima  of  the  skin-friction  magni¬ 
tude  are  used  to  determine  the  separation  locations.  Steady  skin-friction  measurements 
are  obtained  at  fourteen  steady  angles  of  attack.  Unsteady  maneuvers  include  the  ramp 
pitchup  maneuvers  simulated  by  the  DyPPiR.  Mean  wall  static  pressures  are  measured 
at  10  °  and  20  °  angles  of  attack.  Surface  oil  flow  visualizations  are  also  used  in  order  to 
examine  the  steady  surface  flow  topology  and  the  separation  locations  qualitatively.  Both 
steady  and  unsteady  tests  are  performed  for  two  model  configurations.  The  barebody 
(axisymmetric  case)  and  the  body  with  the  sail  (sail-on-side  case). 

Steady  and  the  unsteady  results  obtained  in  this  study  give  valuable  information  about 
the  complex  surface  flow  structure  over  the  DARPA2  submarine  model.  The  steady  skin- 
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friction  values  and  mean  pressure  measurements  are  used  to  describe  the  steady  flow-field 
topology.  Unsteady  data  show  the  difference  between  the  steady  and  the  unsteady  flow- 
fields,  especially  in  the  separation  topology.  In  this  study,  a  first-order  time  lag  model 
is  used  to  approximate  the  unsteady  data.  The  results  of  this  approximation  as  well  as 
the  unsteady  data  may  supply  important  information  to  the  development  of  the  realistic 
flow  models  to  be  used  in  calculating  the  unsteady  3-D  flows  over  complex  geometries. 

The  organization  of  the  chapters  here  can  be  summarized  as  follows.  In  chapter  2, 
the  experimental  apparatus,  facilities  and  techniques  are  described.  This  chapter  gives 
the  details  about  the  DyPPiR,  hot-film  sensors,  constant  temperature  anemometers  and 
the  pressure  measurement  system  used  in  the  experiments.  The  wind  tunnel  model 
including  the  geometry,  different  configurations  and  integration  of  the  sensors  is  also 
described.  Chapter  3  gives  the  test  conditions  and  describes  the  steady  measurements 
and  the  unsteady  maneuvers  performed  in  this  study.  Calibration  of  the  hot-film  sensors 
is  explained  in  chapter  4.  Chapter  5  includes  the  results  obtained  from  the  steady  skin- 
friction  and  pressure  measurements.  The  discussion  about  the  data  and  conclusions  based 
on  the  steady  results  are  presented.  The  steady  data  reduction  procedure  both  for  the  hot- 
fllm  and  the  pressure  measurements  is  also  described  in  this  chapter.  The  unsteady  results 
and  the  data  reduction  procedure  are  presented  in  chapter  6.  The  unsteady  flow  topology 
and  the  differences  between  the  steady  and  the  unsteady  flow  fields  are  discussed.  This 
chapter  also  includes  the  algebraic  and  the  first-order  time  lag  model  approximations  to 
the  unsteady  separation  location  data.  Discussion  about  the  unsteady  data  and  the  time 
lag  approximation  results  are  given.  The  last  chapter  presents  the  overall  conclusions 
obtained  from  this  study.  Uncertainties  in  the  measured  skin-friction  values  and  the 
related  uncertainly  calculations  are  included  in  the  appendices. 


Chapter  1.  Figures 


8 


Figure  1.1:  Limiting  streamline  pattern  and  surfaces  of  separation  for  three  types  of  3-D 
separation:(a)  horseshoe  type  separation;  (b)  Werle  type  separation  (the  view  of  the  surface 
separation  is  rotated  90  °  from  the  view  of  the  limiting  streamline  pattern);  (c)  cross  flow 
separation.  Taken  from  Yates  and  Chapman  [1]. 
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Figure  1.2:  On-axis  cross  sectional  view  of  secondary  flow  streamlines  from  a  crossflow  sep¬ 
aration.  PI:  primary  separation,  P2:  secondary  separation,  Rl:  Primary  reattachment,  R2: 
secondary  reattachment. 


Chapter  2 


Experimental  Facilities,  Equipment 
and  Appciratus 

2.1  Stability  Wind  Tunnel 

The  Virginia  Tech  Stability  Wind  tunnel  is  a  continuous,  closed  jet,  single  return,  sub¬ 
sonic  wind  tunnel  with  six  foot  cross  section  interchangeable  test  sections.  Each  test 
section  is  24  feet  long.  The  tunnel  is  powered  by  a  600  HP  D.C.  motor  driving  a  14  foot 
propeller  [4].  Surface  flow  experiments  on  DARPA2  model  has  been  performed  in  6'  x  6' 
square  test  section.  In  this  conflguration,  a  maximum  speed  of  275  ft/s  and  a  Reynolds 
number  up  to  1.66  x  10®  per  foot  can  be  obtained. 

Tunnel  speed  is  regulated  by  a  custom  designed  Emerson  VIP  ES-6600  SCR  Drive.  This 
drive  system  eliminates  all  the  cyclic  unsteadiness  in  tunnel  velocities  and  turbulence  in¬ 
ducing  vibrations  inherent  with  older  systems.  The  Stability  Tunnel  has  a  low  turbulence 
intensity.  The  small-scale  turbulence  intensity  is  on  the  order  of  0.03%  or  less  depending 
on  the  tunnel  speed.  The  unsteadiness  due  to  large-scale  pulsation  of  the  fan  was  also 
found  to  be  on  the  order  of  0.03%  [4]. 

The  tunnel  temperature  is  stabilized  by  the  use  of  an  air  exchange  tower.  After  a  certain 
amount  of  operatioii  time,  tunnel  temperature  becomes  equal  to  the  outside  temperature. 
Therefore  the  changes  in  the  outside  temperature  directly  effects  the  tunnel  temperature. 
Since  the  hot-film  sensors  are  sensitive  to  the  ambient  flow  temperature  as  well  as  the 
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local  skin-friction  value,  the  changes  in  the  free-stream  temperature  should  be  taken  into 
account  in  the  data  acquisition  and  the  reduction  process. 

Surface  hot-film  measurements  on  DARPA2  model  have  been  performed  in  a  slotted- 
wall  test  section  with  38%  open  air  ratio.  The  main  purpose  for  using  the  slotted-wall 
configuration  in  the  experiments  is  to  reduce  the  blockage  effects  typically  encountered 
during  testing  at  high  angles  of  attack  with  large  models.  More  information  about  the 
slotted  wall  configuration  in  the  stability  tunnel  can  be  found  in  [14].  Figure  2.1  shows 
the  top  view  of  the  Stability  Wind  Tunnel. 


2.2  DyPPiR 

Unsteady  maneuvers  have  been  performed  by  using  the  Dynamic  Plunge,  Pitch  and 
Roll  actuator  of  the  Virginia  Tech  Stability  Wind  Tunnel,  also  known  as  the  DyPPiR. 
The  DyPPiR  was  designed  to  provide  the  unique  capability  of  performing  general,  high 
excursion,  large  scale,  high  Reynolds  number,  unsteady  maneuvers  [4].  The  DyPPiR  has 
three  degrees  of  freedom:  plunge  with  a  range  of  ±0.64  m  measured  from  the  tunnel 
centerline,  pitch  with  a  motion  range  of  ±45  °  and  roll  with  a  range  of  ±140  °.  For  each 
degree  of  freedom,  there  is  an  independent  3000  psi  hydraulic  actuator  which  can  move  a 
maximum  model  load  of  45  kg  and  250  kg  of  hardware  load  at  rates  approaching  9  m/s 
in  plunge  and  120  °/s  in  pitch.  The  DyPPiR  has  also  the  capability  of  performing  ma¬ 
neuvers  around  an  arbitrary  model  center  of  rotation,  which  is  important  for  obtaining 
the  correct  lateral  velocities  across  the  aircraft,  missiles  and  submarines.  The  DyPPiR 
is  computer  controlled  and  in  addition  to  traditional  sinusoidal,  ramp  and  snaking  ma¬ 
neuvers,  any  pre-programmed  real  maneuver  can  be  simulated,  including  time  history 
effects.  Previous  tests  have  confirmed  that  the  sting  mount  places  the  model  far  enough 
from  the  main  DyPPiR  strut  to  make  strut  interference  negligible  [3].  Figure  2.2  shows 
the  DyPPiR  installed  in  the  stability  wind  tunnel.  For  DARPA2  experiments,  besides 
unsteady  maneuvers,  the  DyPPiR  has  also  been  used  as  a  computer  controlled  model 
mount  for  the  steady  skin  friction  measurements.  Due  to  a  technical  problem  in  the 
hydraulically  powered  roll  actuator,  a  manually  controlled  dummy  roll  actuator  has  been 
built  and  used  to  change  the  roll  angle  of  the  model.  This  enabled  a  full  360°  coverage  of 
the  model  roll  orientation.  The  coordinate  nomenclature  for  the  DyPPiR  and  the  model 
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is  shown  in  figure  2.3:  x  is  measured  from  the  nose,  (/>  is  the  circumferential  location 
measured  from  the  windward  line  of  symmetry,  z  is  the  plunge  ordinate,  and  a  is  the 
pitch  angle  and  equivalently  the  model  center  of  rotation  angle  of  attack. 


2.3  Data  Acquisition  System 

The  tunnel  data  acquisition  system  has  been  used  for  sampling  the  analog  signals  from 
each  hot-film  anemometer  as  well  as  recording  tunnel  dynamic  pressure,  static  pressure 
and  the  tunnel  temperature  in  steady  and  unsteady  skin  friction  measurements.  Since 
the  tunnel  data  acquisition  system  is  incorporated  and  synchronized  with  the  DyPPiR 
control  system,  it  has  also  been  used  for  recording  the  plunge,  pitch  and  roll  command 
and  feedback  signals  (voltage  values).  These  values  were  converted  to  proper  dimensional 
units  by  using  a  linear  calibration  routine  which  gave  the  location  and  the  orientation  of 
the  model  as  a  function  of  time. 

The  centerpiece  of  the  system  is  an  SCXI-1001  Mainframe  from  National  Instruments. 
The  mainframe  allows  for  the  installation  of  up  to  12  SCXI  modules  which  may  per¬ 
form  any  desired  signal  conditioning  and  sampling  of  the  input  signals.  Currently  four 
SCXI-1120D  isolation  modules  are  installed  in  the  mainframe.  Each  of  these  modules 
allows  isolation  and  amplification  of  8  differential  analog  voltages  with  an  input  range  of 
-10  to  -1-10  volts.  Each  module  multiplexes  the  8  inputs  which  are  read  by  a  National 
Instruments  AT-MIO-16-XE-10  Data  Acquisition  Card  installed  in  the  tunnel  computer. 
This  results  in  a  total  of  32  isolated  differential  analog  input  channels  with  an  Analog  to 
Digital  conversion  resolution  of  16  bits.  The  software  used  for  data  acquisition  is  written 
using  LabView  4.0  under  the  Windows  NT  environment. 


2.4  Skin-Priction  Measurement  System 

Hot-film  sensors  mounted  on  the  surface  of  the  DARPA2  model  were  used  to  measure 
steady  and  unsteady  skin  friction  magnitudes  in  the  experiments.  Each  hot-film  sensor 
was  connected  to  a  constant  temperature  type  anemometer.  Hot-film  sensors  and  the 
constant  temperature  anemometers  were  the  same  as  the  ones  used  by  Wetzel  [2]  in  his 
prolate  spheroid  study. 
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2.4.1  Hot-film  Sensors 


Theory  of  Operation 


The  thermal  or  hot-film  sensor  benefits  from  the  fact  that  the  heat  transfer  from  a 
sufficiently  small  heated  surface  depends  only  on  the  flow  characteristics  in  the  viscous 
region  of  the  boundary  layer  [15].  Due  to  the  similarity  between  gradient  transport 
of  momentum  and  scalars  (heat),  the  amount  of  heat  transfer  into  the  fluid  gives  a 
measure  of  the  wall  shear  stress  r^.  The  hot-film  gauge  consists  of  a  thin  metallic  film 
positioned  into  a  substrate.  Usually  the  gauge  forms  one  part  of  a  Whetstone  bridge 
in  the  anemometer  circuit  and  an  electric  current  is  passed  through  the  film  in  order  to 
maintain  it  at  constant  temperature  (if  a  constant  temperature  anemometer  is  used)  as 
heat  is  continuously  being  transferred  from  the  film  to  the  moving  fluid  as  well  as  to  the 
film’s  substrate. 


The  ohmic  (joule)  heating  in  the  device  Qj  is  transferred  both  to  the  fluid  and  to  the 
surrounding  substrate.  This  can  be  expressed  as: 

Qj  =  Qs  +  Qf  (2-1) 

where  Qf  represents  the  average  heat  transferred  to  the  fluid  directly  from  the  heated 
surface  and  indirectly  through  the  heated  portion  of  the  substrate.  Qs  represents  the 
average  heat  lost  irretrievably  to  the  substrate. 

The  relationship  between  the  heat  transfer  Q /  and  the  wall  shear  stress  can  be  obtained 
by  making  an  unheated  starting  length  forced  convection  analysis  from  a  single  heated 
sensor  [7].  It  is  assumed  that  the  thermal  boundary  layer  developed  on  the  sensor  is 
within  the  viscous  sublayer  of  a  pre-existing  momentum  boundary  layer. 


We  start  our  derivation  for  the  relationship  between  Qf  and  with  the  thermal  energy 
transport  equation  near  the  wall  for  locally  two  dimensional  flow: 


dT  d^T 
^  dx  ^  dy'^ 


(2.2) 


Here  x  is  the  streamwise  coordinate  and  y  represents  the  normal-to-wall  coordinate.  The 
temperature  is  T,  U  is  the  velocity  in  the  x  direction,  and  a  is  the  thermal  diffusivity. 
Since  the  flow  upstream  of  the  sensor  is  unheated,  the  incoming  flow  temperature  is  the 
same  as  the  free-stream  temperature  Too.  We  assume  the  sensor  length  in  the  streamwise 
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direction  is  short  enough  so  that  the  thermal  boundary  layer  grows  from  the  upstream 
edge  of  the  sensor,  but  remains  entirely  within  the  near  wall  region.  In  the  viscous 
sublayer,  we  can  approximate  U  as: 


U  = 


(2.3) 


where  j  is  the  mean  velocity  gradient  at  the  wall  and  the  shear  stress  at  the  wall 
is  given  by: 


du\ 
M  Hr 


(2.4) 


\9yJ^ 

Here  fi  is  the  dynamic  viscosity  of  the  fluid.  By  putting  non-dimensional  temperature  T 


T  -T 

rp  _ 

^  T  —T 

into  equation  2.2  and  using  equation  2.3,  we  will  have: 

fdU\  df  d'^f 


(2.5) 


(2.6) 


The  boundary  conditions  for  this  problem  in  terms  of  T  would  be: 


r  =  0  at  y  =  0 

(2.7) 

T  =  1  as  y  oo 

(2.8) 

T  =  1  at  X  =  0 

(2.9) 

To  reduce  the  partial  differential  equation  given  by  2.6  to  an 
tion,  we  introduce  the  similarity  parameter  r]  defined  as: 

ordinary  differential  equa- 

^  (3ox)^"^ 

(2.10) 

And  the  resulting  ordinary  differential  equation  would  be: 

f"  -H  rff  =  0 

(2.11) 

Boundary  conditions  for  equation  2.11  in  terms  of  77  become: 

0 

II 

0 

II 

(2.12) 

T  =  1  as  rj  ^  00 

(2.13) 

2.4  Skin-Friction  Measurement  System 


The  solution  of  the  ordinary  differential  equation  can  be  obtained  as: 


T  =  M 


(2.14) 


where 


M-i  =  =  ir  Q  j  =  0.893 


The  heat  flux  at  the  surface  is  given  by: 


Qw  ~~  ^ 


Here  k  is  the  thermal  conductivity.  By  using  equation  2.14  and  the  chain  rule; 


we  obtain  as: 


^ 

dy  drj  dy 


=  Mk{Too  -  Tu,)- 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


oo  2/3 

(Saa;)  ' 

Assuming  that  we  have  a  rectangular  sensor  of  constant  width  w  in  the  spanwise  direction 
and  length  of  I  in  the  streamwise  direction,  we  can  obtain  the  average  heat  transfer  rate 
Qfhy  integrating  qyj  over  w  and  l\ 


I  /  ou 

r  3  \  9i/  > 

Qf  =  w  q-u,dx=  -Mk(T^  -  Too)- — - 
Jo  ^  (3a) 


Then  by  definition  mean  film  coefficient  h  would  be  equal  to: 

h  = 

{Tn.  -  Too)  2  (3a0  ^ 


(2.19) 


(2.20) 


As  can  be  seen  from  equation  2.20,  h  is  directly  proportional  to  >  thus  the  wall 

shear  stress  r^,.  All  the  other  parameters  are  properties  of  the  fluid  and  the  sensor.  T^ 
represents  the  temperature  of  the  hot-film  sensor  and  in  case  a  constant  temperature 
anemometer  circuit  is  used,  this  temperature  remains  constant.  In  the  ideal  case,  the 
free-stream  temperature  Too  can  also  be  considered  as  constant  although  in  most  cases, 
especially  in  long  experimental  run  periods,  there  may  be  considerable  changes  in  the 
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flow  temperature.  (This  issue  and  the  effects  of  free-stream  temperature  change  to  the 
calibration  of  the  hot- film  sensors  will  be  discussed  in  chapter  4). 


If  we  assume  that  the  heat  transfer  to  the  substrate  is  minimized  by  using  isolation 
methods  and  negligible,  then  in  equation  2.1,  the  joule  heating  Qj  will  be  equal  to  the 
heat  transfer  to  the  fluid  Qf.  For  the  hot-film  sensor,  we  can  write  the  joule  heating  term 

p2 

Qj  =  Qj  =  -g-  (2.21) 

■Tlyj 

where  E  is  the  voltage  value  across  the  hot-film  sensor,  and  Rw  is  the  resistance  of  the 
hot-film  sensor  at  its  operating  temperature  T^j,.  Since  is  constant  for  a  constant 
temperature  anemometer,  then  will  also  be  constant  by  considering  the  well-known 
equation: 

R^  =  Rc[l  +  aR  (T^  -  To)]  (2.22) 


In  the  above  equation,  Rc  is  the  resistance  of  the  hot-film  in  a  known  temperature  Tc 
and  an  is  the  temperature  coefficient  of  resistivity.  By  using  equation  2.20  in  accordance 
with  the  equation  2.21; 

(2.23) 


T  —T 

■‘-w  ■‘■oo 


oc  r. 


w 


Therefore,  after  a  proper  calibration,  one  can  determine  the  shear  stress  value  Tu,  by 
measuring  E. 


Specifications  of  the  Hot-film  Sensors  used  in  the  Experiments 

To  measure  the  skin-friction,  hot-film  sensors  designed  and  documented  by  Simpson  et 
al.  [7]  were  used.  Figure  2.4  shows  the  top  view  of  a  typical  sensor.  The  hot-film  sensors 
are  made  of  Balco  foil  (70%  nickel,  30%  iron)  with  a  nominal  temperature  coefficient  of 
resistivity  of  0.0051/  °C  and  are  manufactured  by  the  MINCO,  Inc.  The  main  sensing 
part  of  the  sensor  is  a  spiral  of  5.1  mm  in  diameter  which  is  approximately  0.23%  of 
the  model  length.  At  the  constant  diameter  section  of  the  model,  the  sensor  occupies 
2.184  °  of  the  model  surface  along  the  circumference.  The  foil  sensor  is  bonded  to  a 
kapton  substrate  and  the  total  sensor  thickness  is  0.0635  mm.  Sensors  are  connected  to 
the  constant  temperature  anemometers  by  22  AWG  wires,  which  are  soldered  on  the  two 
tabs  of  the  foil. 
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It  should  be  noted  that  the  hot-film  sensors  have  been  designed  to  measure  only  the 
magnitude  of  the  skin-friction  as  described  in  Simpson  et  al.  [7].  However,  in  his  prolate 
spheroid  study,  Wetzel  [2]  detected  a  directional  sensitivity  of  as  much  as  ±5%.  He 
reduced  the  directional  sensitivity  effects  to  less  than  3%  by  mounting  the  sensors  at 
65  °  off  the  longitudinal  line  on  the  model  surface.  In  this  study,  the  sensor  were  placed 
parallel  to  the  longitudinal  line  on  the  model  surface.  Wetzel  [2]  also  showed  that  the 
sensors  are  incapable  of  resolving  small  scale  turbulence  structures  due  to  the  relatively 
large  sensor  size.  The  skin-friction  magnitudes  measured  by  the  hot-film  sensors  are 
actually  the  values  spatially-averaged  over  the  sensing  spiral  part. 

2.4.2  Constant  Temperature  Anemometers 

The  hot-film  sensors  were  operated  with  Miller-type  non-linearized  constant  temperature 
anemometers  [16].  The  original  boards  were  modified  by  Wetzel  [2]  in  the  prolate  spheroid 
study.  For  the  Suboff  tests,  some  further  adjustments  were  made  in  order  to  increase 
the  stability  of  the  anemometers  and  to  reduce  the  noise  associated  with  the  electronics. 
The  main  power  supply  was  replaced  with  a  new  one  that  can  supply  sufficient  amount 
of  current  to  each  anemometer  at  long  run  periods  and  at  relatively  high  speeds.  This 
change  enabled  the  use  of  fifteen  anemometers  simultaneously  during  the  tests  without 
any  power  loss.  Some  of  the  op-amps  in  the  anemometer  circuits  that  were  not  functioning 
properly  were  replaced  with  the  new  ones  and  this  reduced  the  noise  level  significantly. 
All  the  channels  were  grounded  properly.  Some  of  the  current  knobs  were  replaced  with 
the  new  ones  in  order  to  increase  the  accuracy  of  setting  the  overheat  ratio.  Constant 
temperature  anemometers  had  a  bridge  ratio  of  60  and  the  hot-film  sensors  were  operated 
at  a  nominal,  but  imprecisely  set,  overheat  ratio  of  1.10. 

The  frequency  response  of  the  Miller-type  anemometers  with  the  present  hot-film  sensors 
was  determined  to  be  approximately  200  Hz  by  Wetzel  [2] .  He  used  the  method  described 
in  Wood  [17]  for  finding  the  frequency  response.  While  this  frequency  response  is  too 
low  to  measure  the  fluctuating  turbulence  quantities,  it  is  high  enough  to  resolve  the 
time-history  of  the  spatially-averaged  skin-friction  values  in  the  unsteady  maneuvers. 
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2.5  Wind  Tunnel  Model 

The  light  weight  Darpa2  model  used  in  the  experiments  has  a  generic  Suboff  undersea 
vehicle  geometry  with  a  scaled  length  (L)  of  2.24  m.  The  model  has  a  bow  region 
for  0.0  <  {x/L)  <  0.23,  a  constant  diameter  region  for  0.23  <  [x/L)  <  0.75  and  an 
afterbody  (stern)  region  for  0.75  <  {x/L)  <  1.0.  The  sail  can  be  detached,  resulting  in 
the  axi-symmetric  configuration  of  the  model  which  can  be  thought  as  a  missile  or  to 
a  certain  extent  an  aircraft  fuselage  geometry.  Geometric  equations  defining  the  shape 
of  the  model  body  and  the  sail  are  given  in  a  DTRC  (David  Taylor  Research  Center) 
report  [18].  The  same  equations  are  also  included  in  Whitfield  [13].  Figure  2.7  shows  the 
computer  generated  image  of  the  Suboff  model  by  using  these  equations. 

The  geometry  and  the  structural  components  of  the  model  (except  the  ring  assembly 
integrated  to  the  model  structure  for  moving  the  sail)  are  exactly  the  same  as  that  of 
the  Darpa2  model  used  in  the  force  and  moment  measurements.  Whitfield  [13]  gives  an 
extensive  description  of  the  steps  followed  during  the  construction  of  the  model  and  the 
structural  details.  The  model  is  mainly  made  of  composite  materials.  Between  the  inner 
and  the  outer  carbon  fiber  layer,  the  vinyl  foam  and  the  aluminum  bars  having  hollow 
square  cross  sections  were  placed  on  the  constant  diameter  region  of  the  model.  On  the 
nose  and  the  stern  regions,  the  honeycomb  material  was  used  to  maintain  the  required 
curvature  between  the  inner  and  the  outer  skin.  The  model  weighs  approximately  10  lbs. 
and  has  a  high  strength  to  weight  ratio.  This  is  an  important  property  required  in  the 
unsteady  maneuvers  both  for  the  structural  integrity  of  the  model  and  the  performance 
of  the  DyPPiR. 

The  model  is  mounted  to  the  DyPPiR  with  a  light-weight  composite  sting  which  was 
also  built  for  the  Darpa2  force  and  moment  measurements  [13].  The  last  15  cm  of  the 
stern  region  of  the  model  was  removed  in  order  leave  enough  space  for  the  sting  and  the 
wires  coming  from  the  sensors.  The  overall  weight  of  the  model  with  the  composite  sting 
is  approximately  30  lbs. 

2.5.1  Sensor  Configuration  and  the  Locations 

Figure  2.8  shows  the  model  geometry  with  the  hot-film  sensor  locations.  More  sensors 
have  been  used  near  the  body-sail  junction  region  to  resolve  the  complex  structure  of  the 
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separated  flow.  Oil  flow  visualization  pictures  also  have  been  used  in  order  to  determine 
the  optimum  sensor  locations.  A  total  number  of  29  sensors  were  used  on  the  Darpa2 
Suboff  body.  Since  only  15  anemometers  could  be  run  simultaneously,  the  sensors  were 
separated  into  two  groups  each  having  14  sensors.  These  sensor  sets  were  labeled  as 
Sensor  Set  A  and  Sensor  Set  B.  Two  ribbon  cables,  each  carrying  the  wires  connected 
to  the  sensors  of  each  set,  were  used.  These  ribbon  cables  were  lead  outside  of  the 
model,  to  the  carriage  part  of  the  DyPPiR  where  they  were  connected  to  the  15  BNC 
cables  coming  from  the  anemometers  via  gold  plated  pin-connectors.  BNC  cables  were 
connected  to  a  single  female  connector.  Switching  from  one  set  of  sensors  to  the  other  was 
simply  done  by  connecting  the  pin-connector  of  the  ribbon  cable  to  be  used  to  the  single 
female  connector.  One  sensor  was  individually  connected  directly  to  an  anemometer  and 
used  as  a  control  sensor.  The  data  from  this  sensor  were  taken  during  all  the  runs 
regardless  of  the  sensor  set  used  in  order  to  check  the  consistency  of  the  measurements 
taken  with  different  sensor  groups  in  each  model  position  and  the  orientation.  The 
barebody  measurements  were  performed  by  using  the  sensors  of  Set  B  and  the  individual 
control  sensor  located  at  certain  positions  on  the  long  row  (figure  2.8).  The  same  sensors 
were  used  in  the  sail-on-side  case  throughout  the  whole  circumferential  locations.  The 
remaining  sensors  (set  A)  were  used  near  the  sail  region  in  order  to  resolve  the  flow 
structure.  The  nominal  circumferential  distance  between  the  long  and  the  short  sensor 
row  is  A0  =  11  °.  At  some  sensor  locations,  this  distance  vary  by  ±1  °.  The  uncertainty 
in  the  measurement  of  the  peripheral  distances  is  ±0.44  °.  Table  2.1  gives  the  xjL 
location  of  each  sensor.  Wetzel  [2]  determined  that  to  avoid  the  heating  interference 
between  adjacent  and  upstream  sensors,  they  must  be  spaced  at  least  1.25  inches  apart 
from  each  other.  The  minimum  distance  between  the  sensors  used  in  the  DARPA2 
experiments  is  1.32  inches  which  indicates  the  lack  of  any  heating  interface.  For  the  sail- 
on-side  configuration,  the  chordline  of  the  sail  is  aligned  with  the  circumferential  location 
(j)  =  270  °.  At  its  maximum  thickness  location,  the  sail  extends  from  (f)  =  262  °  to  278  ° 
and  is  placed  between  xjL  =■  0.21  and  x/L  =  0.31. 

2.5.2  Sensor  Mounting  Strategy 

Instead  of  directly  mounting  the  sensors  on  the  model  surface,  the  hot-film  sensors  were 
first  glued  onto  the  cylindrical  plugs  made  up  of  epoxy  and  the  hardener  as  shown  in 
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figure  2.5.  The  upper  surface  of  the  plugs  (where  the  sensors  were  mounted)  were  shaped 
so  that  they  match  with  the  contour  of  model  surface.  The  plugs  have  circular  holes  just 
beneath  of  the  sensing  spiral  part  of  the  sensors.  These  holes  were  filled  with  insulating 
foam  to  minimize  the  amount  of  the  heat  transfer  into  the  model  skin.  Behind  these, 
additional  holes  for  the  pressure  taps  were  drilled  in  the  plugs.  At  the  sensor  locations, 
hollow  circular  plastic  cylinders  were  integrated  to  the  model  structure  (figure  2.9).  The 
sensor  plugs  were  inserted  into  these  holes  (that  had  approximately  the  same  diameters 
as  the  plugs)  and  were  secured  by  using  the  hot-glue.  This  sensor  installation  strategy 
shown  in  figure  2.6  not  only  improves  the  surface  quality  in  the  vicinity  of  the  sensors 
but  also  keeps  the  structural  integrity  of  the  model.  The  plug  region  where  the  sensor 
tabs  were  glued  is  inclined  to  the  inner  part  of  the  model  so  that  the  wires  that  were 
soldered  to  the  sensor  tabs  remain  under  the  surface  level.  The  resulting  gap  was  filled 
with  clay  which  was  shaped  to  be  levelled  with  the  model  surface. 

2.5.3  Rotatable  Ring 

The  sail  was  mounted  on  a  movable  ring  integrated  into  the  model  structure  and  can  be 
rotated  and  fixed  to  any  desired  circumferential  location  with  2°  increments  (figure  2.9). 
The  ring  is  secured  at  the  desired  position  by  a  set  screw.  Both  in  steady  and  unsteady 
testing,  the  skin  friction  magnitude  distribution  on  the  whole  surface  of  the  model  with 
desired  circumferential  spacing  can  be  obtained  by  using  only  one  row  of  sensors  on  the 
model  surface  (note  that  the  purpose  of  the  second  short  sensor  row  shown  in  figure  2.8  is 
to  resolve  the  skin-friction  distribution  around  the  sail  in  the  x  direction).  The  model  is 
rotated  with  a  certain  roll  angle  in  one  direction  while  the  sail  is  rotated  in  the  opposite 
direction  with  the  same  amount.  By  this  action,  while  changing  the  circumferential 
location  of  the  row  of  sensors,  the  model  geometry  and  the  alignment  relative  to  the 
free-stream  for  a  specific  pitch  angle  and  sail  location  is  kept  constant.  A  small  gap 
between  the  model  and  the  sail  has  been  kept  in  order  the  sail  to  be  moved  on  top  of 
the  sensors  without  any  contact.  This  gap  has  been  sealed  properly  by  using  scotch  tape 
and  without  giving  any  fillet  effect  for  every  roll  angle  before  the  data  were  taken. 
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2.5.4  Boundary  Layer  Trips 

In  order  fix  the  transition  location  on  the  model  and  further  guarantee  a  less  Reynolds- 
number  sensitive  separation,  trip  posts  of  0.76  mm  high  cylinders  with  1.28  mm  diameter, 
2.5  mm  spaced  apart,  were  placed  on  the  model  nose  part  at  a;/L  =  0.10.  Two  rows  of 
the  same  trip  posts  were  put  on  the  lower  and  the  upper  surfaces  of  the  sail  in  the 
spanwise  direction.  The  rows  were  located  0.64  cm  (measured  on  the  surface)  away  from 
the  leading  edge  of  the  sail.  The  same  trips  were  used  in  the  prolate  spheroid  study  of 
Wetzel  [2]  and  in  the  DARPA2  force  and  moment  study  of  Whitfield  [13]. 


2.6  Pressure  Measurement  System 

Mean  static  pressures  on  the  model  surface  were  measured  at  29  port  locations  listed  in 
table  2.1.  Each  location  is  behind  a  hot-film  sensor  on  the  same  sensor  plug,  and  the 
distance  between  the  pressure  port  and  the  hot-film  sensor  (measured  from  the  center  of 
the  sensing  spiral  part)  is  approximately  0.2%  of  the  overall  model  length  (figure  2.5). 
Each  pressure  tap  of  each  station  was  connected  to  a  Scanivalve  diaphragm  (having  48 
ports)  through  tygon  tubing.  These  diaphragms  were  connected  to  a  Scanivalve  system 
(CTLR2P/S2-S6  Scanivalve  Corp.)  that  has  a  48  to  1  multiplexing.  In  order  to  reduce 
the  length  of  the  tygon  tubing  that  connects  each  pressure  tap  with  the  scanivalve  system, 
the  housing  of  the  scanivalve  was  placed  inside  the  model  and  mounted  on  the  sting  by 
using  cable  tie-wraps.  The  pressures  were  sensed  by  two  pressure  transducers,  a  Setra 
239  pressure  transducer  with  a  calibrated  range  of  0.0  to  15.0  inches  of  water  was  used  to 
measure  the  free-stream  dynamic  pressure  (To,oo  ~  Poo),  and  another  Setra  239  pressure 
transducer  with  a  calibrated  range  of  -2.5  to  2.5  inches  of  water  was  used  to  measure 
the  static  pressure  relative  to  the  free-stream  static  pressure  {P  -  Poo)  at  the  pressure 
taps.  Olcmen  et  al.  [19]  used  an  inclined  manometer  with  a  resolution  of  0.01  inches  to 
verify  the  calibration  of  the  two  pressure  transducers.  The  total  pressure  Po.oo  and  the 
free-stream  static  pressure  Poo  were  measured  by  using  the  tunnel  Pitot-static  tube.  The 
pressure  coefficient  Cp  at  each  port  is  calculated  by: 

~  r,  _  p 

-‘0,00  -Poo 


(2.24) 
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A  data  acquisition  board  (DT2801)  installed  in  a  PC-AT386  computer  was  used  to  acquire 
data  from  the  pressure  transducers  at  a  sampling  frequency  of  1000  Hz.  The  data  were 
sampled  continuously  for  3.0  seconds  after  a  settling  time  of  10.0  seconds. 

For  the  same  pressure  measurement  system,  the  net  uncertainty  in  Cp  was  reported  as 
±0.018  by  Olcmen  et  al.  [19]. 
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Table  2.1;  Hot-film  sensor  locations  (port  number  is  the  order  of  the  sensors  in  xjL  direction 
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Figure  2.1:  Top  view  of  the  Virginia  Tech  Stability  Wind  Tunnel 


Figure  2.2:  Dynamic  Plunge-Pitch-Roll(DyPPiR)  Model  Mount  installed  in  the 
nel. (taken  from  Wetzel  [2]) 
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Figure  2.5:  Sensor  plug  with  a  hot-film  sensor  mounted 


Figure  2.6:  Schematic  of  the  sensor  mounting  strategy 
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Figure  2.7:  Computer  generated  3-D  view  of  the  Suboff  model  with  the  sail 


detachable  sail 


Figure  2.8:  Top  and  side  view  of  the  model  and  the  hot-film  sensor  locations 
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Figure  2.9:  Barebody  and  sail-on-side  configurations  of  the  model 


Chapter  3 


Steady  Measurements  and  the 
DyPPiR  Maneuvers 


All  the  skin  friction  and  the  pressure  measurements  were  done  with  a  tunnel  speed  of 
42.7±1%  m/s.  For  these  experiments,  the  Reynolds  number  based  on  the  model  length 
Rbl  was  5.5  x  10®.  The  oil  flow  visualizations  were  taken  at  Rei  =  4.5  x  10®  by  Wetzel  [20] . 
Figure  3.1  shows  different  circumferential  regions  of  the  Darpa2  model  where  the  skin- 
friction  and  pressure  measurements  were  performed.  The  angle  range  of  each  region  was 
determined  by  examining  the  oil-flow  visualization  pictures  of  the  sail-on-side  case  taken 
at  10°  and  20°  angles  of  attack.  Since  only  15  anemometers  could  be  used  simultaneously, 
hot-film  measurements  for  each  roll  angle  were  done  by  using  one  set  of  15  sensors.  For 
the  barebody  and  sail-on-side  case,  except  regions  IV  and  V  in  figure  3.1,  both  steady 
and  unsteady  hot-film  data  were  acquired  by  using  sensor  set  B.  In  order  to  resolve  the 
surface  flow  spatial  structure  near  the  sail,  all  the  steady  and  unsteady  measurements 
were  repeated  in  regions  IV  and  V  by  using  sensor  set  A  for  the  sail-on-side  case. 


3.1  Skin-Friction  Measurements 


Measurement  of  the  skin  friction  in  steady  experiments  and  unsteady  maneuvers  has 
been  performed  for  two  model  configurations:  barebody  and  sail-on-side  cases.  In  the 
barebody  case,  both  steady  and  unsteady  skin  friction  data  were  acquired  between  0  =  0° 
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and  180  °  in  the  circunaferential  direction  with  10  °  increments  on  the  windward  side  and 
for  every  2  °  on  the  leeward  side.  For  the  sail-on-side  case,  measurements  were  made 
between  </>  =  0  °  and  360  The  <j)  increment  was  again  2  °  on  the  leeward  side  in  order 
to  locate  the  crossflow  separation  locations  with  low  uncertainty.  On  the  windward  side, 
measurements  were  made  with  10  °  increments  except  the  region  between  (/>  =  270  °  and 
292  °  where  Cj  was  measured  every  2  °  so  as  to  resolve  the  surface  flow  structure  in  the 

vicinity  of  the  sail. 


3.1.1  steady  Hot-film  Measurements 

The  steady  skin-friction  measurements  were  done  using  the  DyPPiR  as  the  model  mount. 
Slotted  walls  were  used  as  the  wind  tunnel  wall  configuration.  Table  3.1  outlines  the 
steady  hot-film  measurements. 


Table  3.1:  Darpa2  steady  skin-friction  measurements 


Body  Configuration 

Pitch  Angles  ( °) 

Sensor  Set 

Region  # 

Roll  Angle  Incr.  ( °) 

barebody 

10  °  &  20  °(w/  solid  walls) 

B 

I 

10.0 

barebody 

10  °  20  °(w/  solid  walls) 

B 

II 

2.0 

barebody 

0.9  to  27.6  °  w/  2  °  incr. 

B 

I 

10.0 

barebody 

0.9  to  27.6  “  w/  2  °  incr. 

B 

II 

2.0 

sail-on-side 

0.9  to  27.6  °  w/  2  °  incr. 

B 

I 

10.0 

sail-on-side 

0.9  to  27.6  °  w/  2  °  incr. 

B 

II 

2.0 

sail-on-side 

0.9  to  27.6  °  w/  2  °  incr. 

B 

III 

2.0 

sail-on-side 

0.9  to  27.6  °  w/  2  °  incr; 

A,  B 

IV 

2.0 

sail-on-side 

0.9  to  27.6  °  w/  2  °  incr. 

A,  B 

V 

2.0 

sail-on-side 

0.9  to  27.6  °  w/  2  °  incr. 

B 

VI 

10.0 

Before  the  actual  tests  with  the  DyPPiR,  the  steady  barebody  Cf  data  were  also  taken 
at  10  °  and  20  °  angles  of  attack  with  the  solid  walls.  The  NACA  Strut  of  the  Stability 
Wind  Tunnel  was  used  as  for  the  model  mount  in  these  experiments.  Steady  hot-film 
measurements  with  the  DyPPiR  were  taken  at  14  angles  of  attack  starting  from  0.9  °. 
The  angle  of  attack  increment  was  about  2  °  and  the  last  angle  covered  was  27.6  .  These 
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angles  are  shown  by  solid  square  symbols  in  figure  3.2.  Besides  obtaining  the  steady 
surface  flow  structure  over  the  model,  the  results  of  the  steady  data  at  these  angles  of 
attack  were  also  used  to  construct  the  quasi-steady  data  to  be  used  in  a  first-order  lag 
model. 


3.1.2  DyPPiR  Maneuvers 

Unsteady  maneuvers  were  performed  by  using  the  DyPPiR.  Slotted  walls  were  used  as 
the  wind  tunnel  wall  configuration.  Table  3.2  summarizes  the  unsteady  skin-fnction 
measurements. 


Table  3.2:  Darpa2  unsteady  skin-friction  measurements 


Body  Configuration 

Maneuver 

Sensor  Set 

Region  # 

Roll  Angle  Incr.  (  °) 

barebody 

0.33  s.  pitch-up 

B 

I 

10.0 

barebody 

0.33  s.  pitch-up 

B 

II 

2.0 

sail-on-side 

0.33  s.  pitch-up 

B 

I 

10.0 

sail-on-side 

0.33  s.  pitch-up 

B 

II 

2.0 

sail-on-side 

0.33  s.  pitch-up 

B 

III 

2.0 

sail-on-side 

0.33  s.  pitch-up 

A,  B 

IV 

2.0 

sail-on-side 

0.33  s.  pitch-up 

A,  B 

V 

2.0 

sail-on-side 

0.33  s.  pitch-up 

B 

VI 

10.0 

Unsteady  results  were  obtained  for  the  pitchup  maneuvers.  The  pitchup  maneuver  per¬ 
formed  for  the  present  work  is  a  simple  linear  ramp  from  1  °  to  27°  in  0.33  seconds.  The 
maneuvers  were  performed  with  a  constant  pitch  rate  of  78  °/s  and  the  model  center 
of  rotation  was  at  xjL  =  0.24.  Figure  3.2  shows  the  DyPPiR  pitch  angle  and  plunge 
location  feedback  for  the  pitchup  maneuvers.  The  DyPPiR  pitch  angle  is  also  the  instan¬ 
taneous  angle  attack  measured  at  the  model  center  of  rotation.  Because  of  the  angular 
motion  of  the  model,  the  instantaneous  local  angle  of  attack  varies  linearly  from  the  nose 
to  the  stern  of  the  model,  with  the  nose  being  at  a  lower  angle  of  attack  than  the  model 
center  of  rotation,  and  the  stern  at  a  higher  angle  of  attack  compared  to  the  model  center 
of  rotation.  The  magnitude  of  the  local  induced  increment  in  angle  of  attack  is  only  a 
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function  of  the  distance  from  the  model  center  of  rotation  since  the  pitch  rate  is  constant. 
For  the  pitch-up  maneuvers  performed,  at  x/L  =  0.0  the  induced  angle  of  attack  incre¬ 
ment  is  approximately  -0.5  while  at  x/L  =  1.0  this  angle  increment  is  approximately 
equal  to  1.5  °.  Although  it  is  possible  to  keep  the  model  center  from  moving  vertically 
by  using  the  plunge  actuator  during  a  maneuver,  it  is  inevitable  that  the  model  center 
of  rotation  translates  downstream  during  a  maneuver.  This  downstream  movement  is 
approximately  0.18  m  and  the  induced  velocity  in  the  downstream  direction  is  around 
1.2%  of  the  free-stream  velocity  which  can  be  considered  as  insignificant.  The  objective 
while  performing  the  pitch-up  maneuver  is  to  get  as  abrupt  a  start  and  stop  as  possible. 
However,  the  DyPPiR  has  a  finite  acceleration  and  deceleration  capability.  Note  that 
the  actual  maneuver  starts  at  t'  =  3.00  in  figure  3.2.  For  each  ({>  orientation  of  the  model, 
the  pitchup  maneuver  was  executed  for  10  times.  In  his  prolate  spheroid  work,  Wetzel  [2] 
has  reported  that  10  repetitions  for  a  given  maneuver  at  each  roll  angle  position  of  the 
model  are  enough  for  ensemble  averages  that  have  a  low  uncertainty  level.  Unsteady  skin 
friction  values  at  each  xjL  measurement  station  has  been  calculated  for  each  repetition 
and  the  final  value  was  obtained  by  ensemble  averaging  the  skin  friction  values  as  will  be 
described  in  the  data  reduction  section  in  more  detailed. 
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Table  3.3:  Darpa2  pressure  measurements 


Body  Configuration 

Pitch  Angles  ( °) 

Region  # 

Roll  Angle  Incr.  (  °) 

barebody 

10  k  20 

I 

10.0 

barebody 

10  k  20 

II 

6.0 

sail-on-side 

10  k  20 

I 

10.0 

sail-on-side 

10  k  20 

II 

6.0 

sail-on-side 

10  &  20  ^ 

III 

6.0 

sail-on-side 

0,  10  k  20 

IV 

2.0 

sail-on-side 

0,  10  k  20 

V 

2.0 

sail-on-side 

10  k  20 

VI 

10.0 
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Steady  mean  surface  static  pressures  on  the  Darpa2  model  were  measured  both  for  the 
barebody  and  the  sail-on-side  cases.  The  NACA  Strut  of  the  Stability  Wind  Tunnel 
was  used  as  for  the  model  mount.  The  information  about  the  pressure  measurements 
are  outlined  in  table  3.3.  For  the  barebody  case,  measurements  were  performed  at  two 
angles  of  attack:  a  =  10  °  and  a  =  20  °.  In  order  to  check  the  symmetry  of  the  flow, 
pressure  data  were  also  taken  on  four  circumferential  locations  of  the  model  {(f)  =  0°, 
90  180  °,  and  270  °)  at  0  °  angle  of  attack.  Barebody  results  were  obtained  both  for 

the  solid  and  slotted  wall  configuration,  whereas  the  sail-on-side  data  were  acquired  with 
the  slotted  walls.  Sail-on-side  pressure  measurements  were  done  at  10  "  and  20  °  angles 
of  attack.  In  regions  IV  and  V,  data  were  also  obtained  for  a  =  0  °. 
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Figure  3.1:  Regions  for  Darpa2  Skin  Friction  and  Pressure  Measurements. 


Figure  3.2:  DyPPiR  plunge  and  pitch  feedback  for  pitchup  maneuvers.  Filled  symbols  show 
the  a  locations  for  the  steady  measurements. 
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Calibration  of  the  Hot-film  Sensors 


The  hot-film  sensors  that  are  used  in  the  actual  tests  have  been  calibrated  for  each  roll 
angle  position  of  the  model  at  zero  angle  of  attack  in  the  Stability  Wind  Tunnel.  During 
the  calibrations,  for  the  sail-on-side  case,  axis-symmetric  flow  around  the  model  was 
maintained  by  removing  the  sail.  The  boundary  layer  velocity  profiles  were  obtained  on 
the  constant  diameter  region  of  the  model  at  0°  angle  of  attack  in  the  Open-throat  tunnel 
of  Virginia  Tech.  Boundary  layer  properties  calculated  from  these  velocity  profiles  were 
used  in  the  calibration  procedure.  The  details  about  the  calibration  process  are  presented 
in  this  chapter. 


4.1  Description  of  the  Experimental  Apparatus 

Boundary  Layer  velocity  profile  measurements  were  performed  in  the  Open-throat  tunnel 
of  the  Virginia  Tech  Aerospace  and  Ocean  Engineering  Department.  The  Open-throat 
wind  tunnel  is  a  wooden,  circular,  return  type  tunnel  with  an  open  test  section.  Detailed 
description  of  the  tunnel  can  be  found  on  the  AOE  Web  page.  In  recent  years,  flow  quality 
of  the  tunnel  has  been  improved  by  installing  a  screen  to  make  the  flow  exit  the  nozzle 
more  uniform  and  by  building  a  vent  to  bleed  the  same  amount  of  air  that  is  entrained 
at  the  test  section  [21].  The  valve  on  the  top  of  the  vent  box  was  kept  open  during  all 
the  runs  in  order  to  make  use  of  this  improvement.  One  disadvantage  of  the  tunnel  is 
the  temperature  control.  It  has  been  seen  that  it  was  difficult  to  keep  the  temperature 
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constant,  however  temperature  change  was  not  high  enough  to  effect  the  velocity  profile 
measurements. 

A  support  assembly,  as  shown  in  figure  4.1,  was  build  in  order  to  place  the  model  in 
the  Open-throat  tunnel.  This  assembly  consists  of  two  wooden  supports  each  connected 
to  the  lower  tunnel  walls  through  two  steel  legs.  The  model  was  placed  on  the  wooden 
supports  in  such  a  way  that  the  whole  part  of  the  constant  diameter  region  could  be 
remained  in  the  test  section.  The  centerline  of  the  model  and  the  support  assembly  in 
the  Y-direction  was  aligned  with  the  centerline  of  the  tunnel  in  the  Y  direction.  However, 
center  of  symmetry  of  the  model  was  shifted  up  from  the  center  of  the  tunnel  with  a 
certain  distance  in  order  to  minimize  the  blockage  effects  and  keep  the  flow  axis-symmetric 
around  the  constant  diameter  part  of  the  model.  Keeping  the  flow  axis-symmetric  is  an 
important  task,  since  all  the  theoretical  calculations  used  to  determine  the  skin-friction 
distribution  are  based  upon  this  assumption.  To  find  the  magnitude  of  the  shift  distance 
in  the  Y-direction,  center  of  the  cross-sectional  area  of  the  model  with  the  wooden  block 
and  the  legs  in  the  Y-Z  plane  has  been  calculated  and  this  center  is  aligned  with  the 
center  of  the  tunnel.  As  a  result,  a  shift  distance  of  3.22  inches  in  the  Y-direction  has 
been  found. 

After  the  model  was  placed  on  the  supports,  velocity  measurements  were  made  with  a 
Pitot  tube  in  eight  stations  to  check  the  symmetry  of  the  flow.  These  eight  stations  were 
taken  on  a  circle  of  radius  30  inches  and  each  station  was  45  °  apart  from  the  next  station 
in  the  circumferantial  direction  as  shown  in  figure  4.2.  (Note  that  streamwise  direction 
X  is  into  the  page  in  figures  4.1  and  4.2).  The  results  are  plotted  in  figure  4.3.  As  can 
be  seen  from  this  plot,  the  velocity  difference  in  each  station  is  in  acceptable  limits  and 
flow  can  be  regarded  as  axis-symmetric., 

Boundary  layer  velocity  profiles  were  obtained  by  using  a  boundary  layer  type  pitot  probe 
with  a  tip  diameter  of  0.035  inches.  The  boundary  layer  probe  was  moved  in  the  normal 
direction  to  the  wall  by  using  an  electrically  driven  traverse  mechanism.  One  revolution 
of  the  barrel  in  the  traverse  mechanism  was  equivalent  to  a  0.025-inch  vertical  move.  The 
probe  was  moved  in  the  streamwise  direction  by  hand  on  an  aluminum  bar  placed  over 
the  test  section  and  secured  in  two  ends  on  the  upper  tunnel  walls.  The  pitot  probe  was 
attached  to  a  water  manometer  with  a  resolution  of  0.02  inches. 
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4.2  Boundary  Layer  Velocity  Profile  Measurements 

Boundary  layer  velocity  profiles  on  the  constant  diameter  region  of  the  model  were  taken 
in  four  different  speeds:  60  ft/s  (18.3  m/s),  80  ft/s  (24.4  m/s),  95  ft/s  (28.9  m/s)  and 
no  ft/s  (33.5  m/s).  Measurements  were  made  in  two  stations  located  on  the  first  sensor 
row  for  each  speed.  As  can  be  seen  from  figure  4.4,  the  first  station  is  at  the  upstream 
location  of  the  constant  diameter  region  {xjL  =  0.25)  and  the  second  station  at  the 
downstream  part  (x/L  =  0.59).  Boundary  layer  properties  obtained  from  the  second 
station  were  used  as  the  initial  value  for  the  theoretical  calculations. 

To  resolve  the  typical  regions  of  the  turbulent  velocity  profiles,  (semi-logarithmic  region, 
wake  region  etc.)  the  pitot  probe  was  moved  with  logarithmic  increments  from  the  wall. 
The  readings  were  assumed  to  be  started  from  the  half  of  the  diameter  of  the  probe, 
0.0175  inches  from  the  wall.  Therefore  the  data  from  the  viscous  sublayer  and  some  part 
of  the  logarithmic  region  couldn’t  be  obtained.  In  fact,  as  will  be  discussed  in  the  next 
section,  the  data  obtained  in  the  vicinity  of  the  wall  weren’t  used  in  the  calculation  of  the 
boundary  layer  properties  because  of  the  possible  errors  originating  from  the  wall-probe 
interference. 

Streamwise  velocity  measurements  were  made  in  29  stations  to  determine  the  boundary 
layer  edge  velocity  Ue  distribution  for  each  speed.  First  and  the  last  stations  are  the  same 
as  the  first  and  the  second  stations  used  in  the  boundary  layer  traverses.  From  station  0 
to  station  27,  the  distance  between  the  stations  are  1  inches.  The  results  obtained  from 
these  measurements  are  given  in  figure  4.8.  As  can  be  seen  from  these  plots,  for  all  the 
speeds,  the  change  in  Ue  in  the  streamwise  direction  is  negligible  and  can  be  regarded  as 
constant.The  temperature  of  the  tunnel  and  the  atmospheric  pressure  have  been  recorded 
regularly  during  all  the  runs  in  order  to  update  the  dynamic  viscosity  /i,  and  the  density 
p  of  the  air  which  were  used  in  the  theoretical  calculations. 

4.2.1  Boundary  Layer  Thickness  Determination 

The  boundary  layer  thickness  8  is  defined  as  the  normal  distance  y  from  the  wall  where 
U  =  0.99  X  Ue-  Here  U  is  the  streamwise  velocity  component  and  Ue  is  the  boundary  layer 
edge  velocity.  The  U  vs.  y  values  have  been  recorded  for  all  the  profile  measurements  and 
the  value  of  the  boundary  layer  edge  has  been  determined  by  making  linear  interpolation 
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between  the  appropriate  values.  Values  of  5  at  the  first  and  the  second  stations  for  each 
speed  are  given  in  Table  4.1. 

Boundary  layer  thickness  determination,  both  in  experimental  and  numerical  studies,  is  a 
difficult  issue,  since  near  the  edge  of  the  turbulent  boundary  layer  velocity  change  is  rather 
small  and  the  uncertainties  associated  with  the  measurement  technique  (e.g.  resolution  of 
the  anemometer  or  reading  errors)  may  cause  significant  deviations  in  the  determination 
of  the  S  values.  Therefore,  the  characterization  of  the  boundary  layer  profiles  with  the 
integral  properties  like  displacement  thickness  5*  and  momentum  thickness  6  would  be 
more  appropriate. 

4.2.2  Calculation  of  the  Boundary  Layer  Properties 

In  order  to  calculate  the  boundary  layer  properties  like  5*  and  6,  the  complete  velocity 
profile  should  be  used.  However,  as  previously  discussed,  the  profiles  obtained  exper¬ 
imentally  didn’t  contain  the  viscous  sublayer  and  some  part  of  the  logarithmic  layer. 
Contributions  to  the  calculation  of  the  integral  properties  from  these  regions  is  not  neg¬ 
ligible,  thus  a  proper  method  to  complete  the  profile  data  should  be  sought.  The  first 
step  is  to  find  the  friction  velocity  Ur  for  that  particular  profile. 

Ur  can  be  obtained  by  using  the  logarithmic  region  in  a  turbulent  profile.  In  the  literature 
one  may  see  many  ways  of  finding  Ur-  All  these  methods  make  use  of  the  logarithmic 
region.  In  this  study  curvature  effects  are  also  considered.  In  White  [22],  the  logarithmic 
region  for  a  turbulent  profile  with  curvature  effect  is  given  by  the  equation; 

U+  =  ^ln{r+)  +  A  (4.1) 

where 

U+  =  —  and  Zn(l-h^)  (4.2) 

Ur  V  \  a/ 

Also  for  the  coefficients,  K=0.41  and  A=5.0  have  been  used  (Coles  and  Hirst  [23]).  At 
r+  =  100,  equation  4.1  gives  U'^  =  16.23.  In  equation  4.2,  a  is  the  radius  of  curvature 
and  for  our  case  in  the  constant  diameter  region  of  the  model  this  value  is  5.25  inches. 
For  finding  Ur,  the  point  in  the  velocity  profile  where  both  U  and  y  values  are  satisfied  is 
sought  for  f/+  =  16.23  and  =  100  by  changing  the  value  of  Ur-  Once  the  point  is  found, 
the  corresponding  value  of  Ur  is  taken  as  the  actual  friction  velocity  for  that  particular 
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profile.  Two  equations  used  in  this  iteration  process  are  derived  from  equations  4.1 
and  4.2  and  given  as: 

y.+=ioo  =  5.25  (ew)  (4.3) 

U  =  16.23C/r  (4.4) 

In  equation  4.3,  C  =  750.0  m~^.  Once  the  friction  velocity  is  obtained,  the  viscous 
sublayer  can  be  approximated  up  to  r"*"  =  10  by  the  equation: 

C/+  =  r+  (4.5) 

To  eliminate  the  bad  points  in  the  vicinity  of  the  wall,  actual  experiment  data  with  r"*" 
value  equal  to  or  greater  than  100  are  used  and  the  region  between  r'^  =  10  and  =  100 
is  approximated  by  the  equation: 

C/+  =  7.5333  (r+)^^®  (4.6) 

In  Kays  and  Crawford  [24],  this  equation  is  given  with  (1/7)  as  the  exponent  and  with  a 
different  constant  as  well.  However  for  the  profiles  obtained  in  this  study,  it  has  been  seen 
that  (1/6)  exponent  form  fits  the  data  better  than  (1/7)  exponent  form  for  all  the  cases. 
The  constant  in  equation  4.6  was  obtained  by  making  equal  to  16.23  for  =  100, 
which  agrees  with  equation  4.1  using  the  Coles  and  Hirst  [23]  constants. 

The  above  procedure  has  been  used  for  obtaining  the  complete  profile  for  each  station 
and  speed.  The  results  are  given  in  figures  4.5  to  4.7.  The  logarithmic  region  gets  larger 
as  the  speed  increases  and  particularly  for  a  given  speed  the  logarithmic  region  of  the 
profile  obtained  in  the  second  station  is  again  larger  with  respect  to  the  one  obtained  in 
the  first  station. 

After  the  complete  profile  data  for  each  case  were  obtained,  boundary  layer  properties  5* , 
0  and  the  shape  factor  H  have  been  calculated.  These  properties  and  Reg,  the  Reynolds 
Number  based  on  6  for  each  case  are  presented  in  table  4.1. 

Although  all  the  calibration  steps  described  to  this  point  have  been  applied  to  74  =  60 
ft/s  case  besides  80,  95  and  110  ft/s,  boundary  layer  properties  obtained  at  this  particular 
speed  were  not  used  in  the  remaining  calibration  procedure.  The  velocity  profiles  at  60 
ft/s  showed  transitional  character.  Since  all  the  theoretical  calculations  in  the  calibration 
procedure  were  based  on  fully  turbulent  axis-symmetric  flow  assumption  over  the  model, 
the  results  of  the  60  ft/s  case  were  not  used  in  the  skin-friction  determination.  Therefore 
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Table  4.1:  Boundary  layer  properties  at  Ue  =  80,  95  and  110  ft/s. 


Ue,  ft/s 

Station  # 

U^,  ft/s 

(5,  inches 

5*15 

9/5 

H 

Ree 

80.0 

1 

1.11 

0.26 

0.17 

0.12 

1.44 

1114 

80.0 

2 

0.98 

0.96 

0.15 

0.11 

1.35 

3840 

95.0 

1 

1.33 

0.22 

0.17 

0.12 

1.43 

1127 

95.0 

2 

1.13 

0.96 

0.15 

0.11 

1.34 

4752 

110.0 

1 

1.44 

0.25 

0.16 

0.11 

1.44 

1470 

no  .0 

2 

1.30 

0.85 

0.15 

0.11 

1.36 

4920 

actual  experimental  data  used  in  the  calibration  procedure  consist  of  boundary  layer 
velocity  profile  measurements  acquired  at  80,  95  and  110  ft/s.  As  can  be  seen  from 
table  4.1,  experimental  calibration  C/  values  were  obtained  for  1010  <  Reg  <  4900. 
The  extension  of  the  calibration  range  to  140  ft/s  was  achieved  by  using  a  theoretical 
approach  described  in  section  3.3. 


For  80,  95  and  110  ft/s,  to  determine  the  9  distribution  between  two  measurement  sta¬ 
tions,  the  following  momentum  integral  equation  has  been  used: 

I  -0.2857 


0.03138 


,  9 

R^a  In  I  1  -l-  9.337 — 


dx 


(4.7) 


This  equation  has  been  obtained  by  using  the  approach  in  Kays  and  Crawford  [24]  and 
making  necessary  modifications  to  include  the  transverse  curvature  effects.  Equation  4.7 
which  represents  an  initial  value  problem,  has  been  solved  numerically  by  using  Modified 
Euler’s  Method  in  order  to  determine  the  momentum  thickness  distribution.  As  the 
initial  value  for  the  momentum  thickness,  the  9  value  measured  at  the  second  station  has 
been  used. 


4.3  Calculation  of  the  Skin-Friction  Values 


After  obtaining  the  9  distribution  at  80,  95  and,  110  ft/s,  the  skin  friction  coefficient  at 
each  sensor  location  has  been  calculated  by  using  the  Ludwieg-Tillmann  equation  [25]: 


V 


^  =  0.123  X 
2 


-0.268 


(4.8) 
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Steady  and  unsteady  experiments  over  the  DARPA2  model  have  been  performed  at  a 
nominal  wind  tunnel  speed  of  140  ft/s.  The  wall  shear  stress  at  this  speed,  especially  on 
the  windward  side  of  the  model  at  an  angle  of  attack,  would  be  significantly  high  com¬ 
pared  to  the  values  obtained  at  the  experimental  calibration  speeds.  Since  the  highest 
calibration  speed  that  can  be  reached  at  the  open  throat  tunnel  was  110  ft/s,  the  cal¬ 
culation  of  the  actual  shear  stress  in  the  experiments  would  require  to  approximate  the 
values  well  above  the  skin-friction  magnitudes  obtained  from  the  experimental  velocity 
profile  measurements.  In  order  to  reduce  the  error  in  the  extrapolation,  the  skin-friction 
distribution  at  the  constant  diameter  region  of  the  model  at  140  ft/s  was  approximated 
by  using  the  equation  given  in  White  [26]: 


Cf  =  0.0015  + 


0.20  +  0.016 


(4.9) 


In  the  same  work.  White  [26]  showed  that  equation  4.9  had  an  rms  error  of  ±9%  compared 
with  the  available  data,  which  was  the  lowest  of  any  theory  known  to  the  author  at  that 
time. 


4.4  Calculation  of  the  Calibration  Coefficients 


In  order  to  relate  the  voltage  values  acquired  from  the  constant  temperature  anemometers 
with  the  wall  shear  values,  the  surface  hot-film  version  of  King’s  Law  has  been  used 
(Bruun  [27]): 

—A— =  4  +  (4.10) 

Here,  E  is  the  time-averaged  voltage  value  obtained  from  a  surface  hot-film  sensor  con¬ 
nected  to  a  constant  temperature  anemometer,  stands  for  the  sensor  temperature  and 
Too  for  free-stream  temperature  of  the  flow  in  the  tunnel.  The  purpose  of  the  calibration 
is  to  determine  the  coefiicients  A  and  B  in  equation  4.10.  For  finding  these  coefficients, 
E  and  the  corresponding  r^,  value  obtained  at  80,  95,  110  and  140  ft/s  have  been  used 
to  make  a  linear  regression. 

In  equation  4.10,  the  change  in  the  free-stream  temperature  will  also  cause  a  change  in 
the  calibration  coefficients  A  and  B.  Since  the  temperature  of  the  stability  tunnel  is 
ambient  and  cannot  be  controlled,  the  calibration  procedure  was  repeated  as  the  tunnel 
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temperature  changed.  In  order  to  minimize  the  uncertainty  in  the  skin  friction  measure¬ 
ments  due  to  the  free-stream  temperature  change  for  the  barebody  case,  the  calibration 
coefficients  for  each  sensor  have  been  re-calculated  for  every  roll  angle  position  of  the 
model  before  taking  steady  data  and  performing  unsteady  maneuvers  for  that  specific 
roll  angle.  At  each  roll  position  of  the  model,  the  voltage  values  E  from  each  sensor 
were  acquired  at  0°  angle  of  attack  for  the  speeds  80,  95,  110  and  140  ft/s.  Since  the 
value  corresponding  to  each  speed  at  all  the  sensor  locations  were  known,  the  calibration 
coefficients  A  and  B  for  each  sensor  and  roll  angle  position  could  be  obtained  by  using 
equation  4.10.  The  free-stream  temperature  change  was  at  most  ±0.5  C  °  between  the 
start  of  the  steady  measurements  and  the  end  of  the  unsteady  maneuvers  for  each  roll 
angle  and  this  was  included  in  the  overall  uncertainty  calculations.  For  the  measurements 
with  the  sail,  calibration  coefficients  for  each  sensor  were  re-calculated  approximately  in 
every  10  roll  angle  by  simply  detaching  the  sail  from  the  body  and  applying  the  calibra¬ 
tion  procedure  to  the  barebody.  Between  each  calibration  runs,  the  same  values  of  A  and 
B  coefficients  were  used  for  each  sensor. 


Figure  4.9  shows  the  skin-friction  vs.  xjL  distribution  obtained  as  the  result  of  the 
calibration  procedure  at  the  actual  experiments  in  the  Stability  Wind  Tunnel  at  0°  angle 
of  attack.  The  Reynolds  number  based  on  the  model  length  Rei  is  5.5  x  10®.  In  the 
same  figure,  the  data  are  compared  with  the  Cf  values  obtained  at  David  Taylor  Model 
Basin  (DTMB)  for  Rcl  =  1.2  x  10^  (Huang  et  al.  [28]).  Skin-friction  measurements  at 
DTMB  were  based  on  the  principle  that  shear  stress  on  a  body  in  a  flow  can  be  measured 
with  using  small  obstacles  that  stagnate  the  velocity  field  near  the  surface  to  produce 
a  pressure  rise  that  is  approximately  proportional  to  the  shear  stress.  Huang  et  al.  [28] 
placed  small  blocks  with  certain  dimensions  on  the  model  surface  near  to  the  pressure 
tap  locations.  The  pressure  was  measured  at  the  tap  with  and  without  these  obstacles 
and  the  measurement  differences  were  used  to  compute  the  shear  stress  at  that  location. 
They  reported  an  uncertainy  of  ±0.0002  in  their  skin-friction  measurements.  Skin  friction 
Cf  vs.  x/L  distributions  of  Virginia  Tech  and  DTMB  follow  the  same  trend.  However, 
the  magnitudes  of  the  Cj  obtained  at  Virginia  Tech  are  bigger  than  that  of  DTMB.  This 
is  likely  due  to  the  difference  between  the  Reynolds  numbers  and  the  uncertainty  of  the 
both  measurements. 
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Figure  4.1:  Cross-sectional  view  of  the  model  and  support  assembly  in  the  Y-Z  plane 


Figure  4.2:  Measurement  stations  for  the  flow  symmetry  check 
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Station  # 


Figure  4.3;  Results  of  the  velocity  measurements  at  the  flow  symmetry  check  points 


Figure  4.4:  Planform  view  of  the  sensor  configuration,  boundary  layer  traverse  and  stremwise 
velocity  measurement  stations  (all  dimensions  are  measured  in  inches  and  the  figure  is  not 
drawn  to  scale). 
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Figure  4.5:  Boundary  layer  velocity  profiles  for  Ue  =  80  ft/s. 


Figure  4.6:  Boundary  layer  velocity  profiles  for  [/g  —  95  ft/: 
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r+ 


Figure  4.7:  Boundary  layer  velocity  profiles  for  Ue  =  110  ft/s. 


Figure  4.8:  Streamwise  Velocity  Distributions  for  ail  calibration  speeds. 
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Figure  4.9;  C/  vs.  x/L  at  0°  angle  of  attack  (For  the  measurements  at  Virginia  Tech  Rei  = 
5.5  X  10®,  for  the  measurements  at  David  Taylor  Model  Basin  (DTMB)  Rei  =  1.2  x  10^) 


O  cf_Virginia  Tech 

♦  cf_DTMB 

i 

D  CQZ2 

0 

♦ 

◄ 

^  s 
►  \ 

►  ^  < 

♦o 

V 

Chapter  5 


Results  and  Discussion  for  Steady 
Measurements 

5.1  Steady  Skin-R-iction  Measurements 

5.1.1  Data  Acquisition  and  Reduction 

A  Labview  software  program  was  developed  for  the  acquisition  of  the  steady  skin-friction 
data.  This  software  program  was  used  to  obtain  the  raw  voltage  values  from  the  constant 
temperature  anemometers,  each  connected  to  a  hot-film  sensor  mounted  on  the  model 
surface.  In  order  get  a  single  steady  skin-friction  Cf  value,  10  seconds  of  voltage  data 
were  acquired  at  a  sampling  rate  of  500  Hz.  This  gave  a  total  number  of  5000  samples 
per  one  skin-friction  value. 

The  steady  values  that  were  measured  or  calculated  can  be  defined  as  arrays  of  indexed 
variables  for  the  description  of  the  steady  data  reduction  process.  For  instance,  the 
voltage  value  acquired  from  each  anemometer  can  be  written  as. 

E  =  E{aj,Sk,<f>m)  (^•^) 

j  =  l,...,14  fc  =  l,...,15  and  m  =  I,  ...,nroll  (5.2) 

where  cxj  corresponds  to  each  steady  angle  of  attack;  Sk  to  the  sensor  number  or  location, 
and  4>rn  to  the  roll  angle  position  of  the  model.  The  total  number  of  roll  angles  nroll  is 
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56  for  the  barebody  case  and  117  for  the  sail-on-side  case. 

The  steps  followed  during  the  steady  skin-friction  data  reduction  process  can  be  described 
in  the  items  below: 


1.  For  each  roll  position  and  angle  of  attack,  the  voltage  data  of  the  free-stream  tem¬ 
perature  T(oij,4>m) ,  the  dynamic  pressure  q{oij,(f)m)>  &nd  the  atmospheric  pressure 
p{aj,<f>m)  converted  to  actual  values  by  using  the  calibration  factors. 


2.  The  wall  shear  stress  T^{aj,Sk,  4>m)  is  calculated  by  using  King’s  Law  (equation  4.10) 
for  each  sensor,  at  each  angle  of  attack  and  model  roll  position: 


r"^(Q:j,  Ski  ^Tn} 


A{ski(t>7n) 
B(^Sk  1  (f^m) 


(5.3) 


A{ski  4>m)  and  B{sk,(j>m)  are  the  calibration  coefficients  obtained  for  each  sensor  at 
each  roll  angle  position  of  the  model.  The  skin-friction  coefficient  Cf{aj,Sk,(t>m) 
can  be  calculated  simply  by  dividing  the  wall  shear  stress  by  the  dynamic  pressure: 


Cfi^CXj ,  Ski  0m) 


ry,{ajiSki(f>m) 


(5.4) 


3.  Before  the  determination  of  the  separation  locations,  Cf  =  Cf{aj,Ski(f>m)  vs.  0  = 
<pm  distributions  for  each  sensor  are  obtained  at  each  angle  of  attack.  In  order  to 
filter  out  some  deviant  points  in  the  Cf  vs.  (p  distributions  of  the  sail-on-side  case, 
a  robust  locally  weighted  regression  technique  called  loess  method  (Cleveland  [29]) 
is  used.  This  smoothing  technique  uses  pre-determined  windows  of  0  values  in  the 
regression  process  and  gives  higher  weights  to  the  points  near  to  the  location  where 
the  smoothed  Cj  value  is  sought.  Outliers  in  the  data  set  are  detected  during  the 
fitting  process  and  not  used  in  the  regression.  These  features  enable  to  smooth 
the  data  without  changing  the  original  pattern.  A  mathematical  description  of  the 
loess  method  and  important  parameters  to  be  considered  in  the  smoothing  process 
is  given  in  Appendix  A.  For  further  details  about  the  method,  see  Cleveland  [30]. 

4.  The  separation  locations  at  each  sensor  location  are  determined  from  the  graphs  of 
circumferential  skin-friction  distributions.  The  local  minima  in  Cf  is  used  to  locate 
the  separation. 
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5.1.2  Results  and  Discussion  for  the  Barebody  Case 

For  the  range  of  test  conditions  mentioned  in  chapter  4,  steady  results  over  the  barebody 
show  typical  characteristics  of  the  crossflow  separation.  Figure  5.1  shows  the  oil  flow 
visualization  of  the  near-wall  fluid  over  the  constant  diameter  region  of  the  model  at 
a  =  20  °  in  the  barebody  case  (Wetzel  [20])  and  qualitatively  describes  the  crossflow 
separation  topology.  Skin-friction  lines  converge  along  the  separation  lines  and  from  this 
figure  two  separation  lines  can  be  identified  on  the  leeward  side  of  the  model.  The  one 
closer  to  the  windward  side  is  defined  as  the  primary  separation  and  the  other  as  the 
secondary  separation  line. 

Figures  5.5  to  5.18  give  the  skin-friction  Cf  distribution  vs.  circumferential  location  on 
the  model  surface  0  for  different  x/L  stations  each  at  a  steady  angle  of  attack  a.  They 
show  the  change  of  the  circumferential  Cf  distribution  as  a  function  of  x/L  at  each  a 
value.  Figures  5.19  to  5.32  give  the  Cf  distribution  vs.  ^  for  all  steady  angles  of  attack  at 
each  sensor  location.  From  these  figures,  the  change  of  circumferential  Cf  distribution  as 
a  function  of  a  at  each  x/L  location  can  be  observed.  Barebody  results  are  represented 
by  the  filled  symbols.  Iii  these  figures,  primary  separation  locations  can  be  detected  as 
the  most  windward  minima  of  the  Cf  distribution  and  the  secondary  separation  location 
can  be  obtained  by  finding  the  second  minimum.  Figures  5.33  and  5.34  show  the  primary 
separation  locations  as  a  function  of  a  for  each  sensor  location.  The  primary  separation 
locations  can  be  first  detected  at  the  last  five  stations  x/L  =  0.638,  0.706,  0.774,  0.819, 
and  0.863  for  a  =  3.1  °.  However,  the  primary  separation  locations  can  be  seen  more 
clearly  at  angles  of  attack  starting  from  a  =  5.1  °.  From  figures  5.33  and  5.34,  it  can  be 
seen  that  as  x/L  increases,  the  primary  separation  location  moves  towards  the  windward 
side  at  a  specific  angle  of  attack.  The  same  behavior  can  also  be  seen  as  the  angle  of 
attack  increases;  the  primary  separation  location  moves  towards  the  windward  side  at 
a  specific  x/L  location.  At  x/L  =  0.266,  <^  =  150  °  can  be  located  as  the  first  primary 
separation  at  a  =  11.3 For  this  station,  the  primary  separation  location  moves  to 
(/)  =  134  °  at  a  =  27.6  °.  At  x/L  =  0.819,  the  primary  separation  location  starts  from 
(^  =  147  °  at  a  =  3.1  °  and  moves  to  </>  =  111  °  at  a  =  27.6  °.  Figures  5.35  and  5.36 
show  the  secondary  separation  locations  vs.  a  for  different  x/L  stations.  The  secondary 
separation  can  be  first  detected  at  a  =  11.3  °.  For  this  angle  of  attack,  the  secondary 
separation  location  is  <^  =  141  °  at  x/L  =  0.774  and  (j)  =  143  °  at  x/L  =  0.819.  In 
general,  at  a  certain  angle  of  attack,  the  secondary  separation  location  moves  leeward  as 
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xfL  increases.  Although  no  data  between  the  measurement  stations  are  available,  the 
separation  location  at  each  xjL  location  can  be  used  to  interpolate  the  values  in-between 
and  this  may  give  an  estimate  of  the  primary  and  secondary  separation  lines  on  the  model 
surface. 

From  the  circumferential  skin-friction  distributions,  a  nearly  flat  profile  can  be  seen 
(figure  5.31)  at  the  vicinity  of  the  minimum  point  for  the  last  station,  xjL  =  0.863 
located  on  the  stern  region.  The  flat  profile  extends  from  (j)  =  100  °  to  (j)  =  124  at 
a  =  15.3  This  profile  indicates  separated  low  speed  fluid  in  this  complex  flow  region 
that  makes  the  identification  of  the  true  minimum  difficult.  The  result  is  consistent  with 
the  oil-flow  visualization  pictures  of  this  region.  Figure  5.2  shows  the  low  speed  fluid 
pattern  on  the  stern  region  at  a  =  15  °. 

Figure  5.19  shows  a  bump  in  the  Cf  distribution  on  the  leeward  side  of  the  model  at  the 
first  station  xjL  =  0.110  starting  from  a  =  15.3  °.  In  this  region,  at  a  =  27.6  °,  the 
increase  in  Cf  starts  approximately  from  ^  =  124  °  and  reaches  the  peak  value  at  around 
0  =  136  °  and  returns  back  to  regular  trend  at  =  144  °.  A  similar  pattern  can  be 
observed  also  for  the  stations  x/L  =  0.131  and  0.170.  These  three  stations  are  located 
on  the  nose  region  of  the  model.  Therefore,  this  Cf  pattern  may  imply  a  weak  separation 
and  reattachment  of  the  flow  on  the  nose  region  of  the  model  at  relatively  high  angles  of 
attack. 

Figures  5.37  and  5.38  show  the  comparison  of  the  barebody  primary  separation  locations 
obtained  with  the  hot-film  sensors  and  the  oil-flow  visualizations.  The  angle  of  attack 
is  approximately  15  °  for  figure  5.37  and  20  °  for  the  other.  In  both  figures,  the  oil-flow 
predicts  the  separations  locations  more  windward  compared  to  the  hot-film  results  with 
a  difference  of  approximately  20°.  As  described  in  Wetzel  et  al.  [32],  the  errors  associated 
with  the  separation  locations  of  the  oil-flow  visualization  can  be  attributed  to  the  gravity 
effects  or  direct  interactions  between  the  flowfield  and  the  oil  mixture,  which  tends  to 
pool  near  separations.  As  in  this  study,  these  errors  can  be  very  significant,  therefore 
oil-flows  should  be  used  for  qualitative  interpretations  of  the  flow  only. 

Steady  barebody  Cf  data  obtained  with  the  solid  walls  are  compared  with  the  slotted 
wall  results  in  figures  5.39  through  5.44  at  three  x/L  stations:  0.345,  0.570,  and  0.819. 
The  results  are  presented  for  two  angles  of  attack:  for  the  solid  wall  case,  a  =  10  °  and 
a  =  20  °  and  for  the  slotted  wall  case  a  =  9.3  °  and  a  =  19.4  °.  For  all  the  stations. 
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the  difference  between  the  solid  wall  and  the  slotted  wall  data  can  be  observed  on  the 
windward  side  and  on  the  last  part  of  leeward  side  starting  approximately  from  (j)  =  150°. 
On  these  regions,  the  skin-friction  values  of  the  slotted  wall  case  are  higher  than  the  ones 
obtained  in  the  solid  wall  case.  This  difference  is  more  significant  at  a  =  20  due  to 
the  increased  blockage  effect  of  the  model.  At  both  angles  of  attack,  the  difference  m 
Cf  decreases  as  x/L  increases.  In  other  words,  the  blockage  effect  m  Cj  values  is  more 
obvious  on  the  upstream  constant  diameter  part  of  the  model  as  compared  to  the  stern 
region.  The  primary  and  the  secondary  separation  locations  are  the  same  m  both  wall 
configurations,  thus  the  separation  locations  are  less  sensitive  to  the  blockage  effect. 

5.1.3  Results  and  Discussion  for  the  Sail-on-side  Case 

In  figures  5.5  to  5.32,  open  symbols  stand  for  the  Cf  vs.  distributions  of  the  sail-on- 
side  case.  It  is  more  appropriate  to  evaluate  the  results  of  the  sail-on-side  case  in  two 
separate  regions:  the  region  with  no  sail  (between  (^  =  0°  and  180°)  and  the  region  with 
the  sail  (between  ^  =  180  °  and  360  °).  In  the  first  region,  the  origin  and  the  variation 
of  the  primary  and  the  secondary  separation  lines  as  a  function  of  x/L  and  a  show  the 
same  characteristics  as  defined  for  the  barebody  case.  This  implies  that  the  main  flow 
feature  on  the  non-sail  region  is  the  cross  flow  separation.  Figures  5.45  and  5.46  give 
the  primary  separation  locations  vs.  a  for  each  x/ L  stations  on  the  non-sail  side.  The 
same  distribution  for  the  secondary  separation  locations  are  presented  m  figures  5.47 
and  5.48.  Although  the  general  trend  is  the  same,  the  locations  of  the  primary  and  the 
secondary  separations  are  slightly  different  from  the  barebody  separation  locations.  They 
are  more  leeward  compared  to  the  barebody  locations.  This  difference  can  be  observed  in 
figures  5.49  and  5.50.  Figure  5.49  shows  the  variation  of  the  primary  separation  location 
with  a  for  three  different  x/L  stations.  As  can  be  seen  from  this  figure,  the  primary 
separation  locations  of  the  barebody  case  and  the  non-sail  region  of  the  sail-on-side  case 
are  approximately  the  same  within  uncertainties  up  to  a  =  15.3°.  Beyond  this  angle,  the 
separation  locations  of  the  sail-on-side  case  start  to  deviate  from  the  barebody  results 
having  an  offset  in  the  leeward  direction.  This  difference  is  obvious  for  a  =  21.4  °,  23.4  ° 
and  25.5°.  Figure  5.50  shows  the  variation  of  the  secondary  separation  location  with  a  for 
the  same  x/L  locations.  In  this  figure,  for  all  angles  of  attack,  the  secondary  separation 
locations  on  the  non-sail  region  of  the  sail-on-side  case  are  shifted  in  the  leeward  direction 
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compared  to  the  secondary  separation  locations  of  the  barebody  case. 

The  skin-friction  distribution  trend  of  the  sail  region  on  the  nose  part  of  model  is  approx¬ 
imately  the  same  as  the  one  obtained  for  the  barebody  case.  At  stations  xjL  —  0.110, 
0.131  and  0.170,  the  same  flow  separation  and  reattachment  pattern  can  be  seen  starting 
from  a  =  15.3  °.  As  can  be  seen  from  the  figures  5.19,  5.20,  and  5.21,  (7/  vs.  0  distri¬ 
butions  for  all  angles  of  attack  are  symmetric  with  respect  to  0  =  180  The  low  speed 
separated  flow  region  on  the  stern  of  the  model  can  be  seen  again  from  figure  5.4  which 
shows  the  oil  flow  pattern  of  this  region  for  the  sail-on-side  case  at  a  =  10 

Downstream  of  the  sail,  the  flow  structure  on  the  sail  side  of  the  model  is  much  different 
than  the  one  observed  for  the  non-sail  side.  The  flow  field  in  this  region  is  strongly  affected 
by  the  presence  of  the  sail.  Compared  to  the  separation  topology  of  the  barebody  and 
the  non-sail  region  of  the  sail-on-side  case,  the  separation  location  trend  as  a  function 
of  x/L  and  a  shows  significant  differences.  Two  minima  in  Cf  vs.  (j)  distributions  on 
the  leeward  side  of  the  sail  region  can  be  observed  for  certain  angles  of  attack  and  x/ L 
locations. 

Figures  5.51  and  5.52  show  the  first  (measured  from  0  =  180  °)  of  these  separation 
locations  as  a  function  of  a  for  different  streamwise  measurement  stations  starting  from 
x/L  =  0.434.  At  the  stations  upstream  of  x/L  =  0.434,  this  separation  pattern  is  not 
observed.  For  the  other  stations,  first  measured  separation  location  on  the  leeward  side 
does  not  follow  a  certain  pattern  as  the  angle  of  attack  or  x/L  are  increased.  It  changes 
between  (f)  =  207  °  and  </>  =  216  °  in  an  irregular  fashion.  This  separation  pattern  can 
be  observed  at  different  angles  of  attack  for  different  sensor  locations.  The  first  angle  of 
attack  at  which  this  separation  can  be  detected  is  a:  =  5.1  °  for  all  sensor  locations,  while 
the  highest  angle  of  attack  changes  for  different  x/L.  Among  all  the  sensor  locations, 
the  highest  angle  of  attack  is  17.4  °  and  observed  at  x/L  =  0.434.  Beyond  this  angle  of 
attack,  the  minimum  associated  with  the  first  separation  location  on  the  leeward  side  of 
the  sail  region  vanishes  at  all  x/L  locations. 

A  second  minima  in  Cf  distribution  can  be  located  on  the  leeward  side  of  the  sail  region. 
Figures  5.53  and  5.54  show  these  separation  locations  at  angles  of  attack  between  a  =  3.1° 
and  a  =  27.6  °  for  the  sensor  stations  starting  from  x/L  =  0.434.  For  the  stations 
x/L  =  0.434,  0.501,  0.570,  and  0.638,  the  separation  locations  move  leeward  as  the  angle 
of  attack  is  increased  up  to  15.3  °.  After  this  particular  angle  of  attack,  the  separation 
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location  for  each  sensor  remains  approximately  the  same.  For  the  stations  downstream 
of  x/L  =  0.638,  the  separation  locations  move  leeward  until  cv  =  9.3  °  is  reached.  At  the 
remaining  angles  of  attack,  the  separation  location  remains  approximately  at  the  same 
(j)  location  for  each  sensor  station.  The  results  on  the  sail  side  indicate  that  the  flow 
field  does  differ  from  the  crossflow  separation  structure  observed  for  the  barebody  and 
non-sail  region  of  the  sail-on-side  case.  Therefore  the  categorization  of  the  separation 
locations  as  the  primary  or  the  secondary  is  not  clear  and  may  not  reflect  the  real  flow 
structure  of  this  region. 

The  flow  in  the  vicinity  of  the  sail-body  junction  is  dominated  by  the  horseshoe  type 
separation.  This  can  be  clearly  seen  from  figure  5.3  which  shows  the  oil-flow  pattern  in 
the  leeward  side  of  the  sail  region  at  a  =  15  °.  The  separation  line  emanating  from  a 
three-dimensional  stagnation  point  upstream  of  the  sail  extends  from  both  leeward  and 
windward  side  of  the  sail  and  travels  downstream.  The  separation  line  on  the  leeward 
side  of  the  sail  can  be  seen  along  the  converging  skin  friction  lines  in  figure  5.3.  Although 
the  flow  topology  near  the  onset  of  the  horseshoe  separation  is  different  from  that  of 
the  crossflow  separation  as  described  in  Yates  and  Chapman  [1],  both  show  similar  char¬ 
acteristics  of  separated  flows  downstream;  there  is  a  strong  convergence  of  the  limiting 
streamlines  on  the  surface  and  there  are  concentrated  regions  of  vorticity  in  the  flow. 
This  may  raise  the  ambiguity  about  the  identification  of  the  separation  pattern  down¬ 
stream  of  the  sail.  However,  the  first  and  second  separation  locations  on  the  leeward  side 
of  the  sail  region  described  in  the  previous  paragraph  may  represent  the  separation  lines 
emanating  fi:om  the  leeward  and  windward  side  of  the  sail  respectively.  At  high  angles 
of  attack,  only  one  separation  line  (described  by  the  second  minimum  in  the  previous 
paragraph)  can  be  detected.  This  may  imply  that  leeward  and  windward  separation  lines 
emerge  just  at  the  downstream  of  the  sail  forming  a  single  separation  line. 

Figures  5.55  to  5.68,  show  the  C/  contours  around  the  sail  region  at  steady  angles  of 
attack.  Note  that  the  blank  areas  are  the  regions  where  no  data  were  acquired.  The 
vertical  blank  between  x/L  ~  0.22  and  x/L  ~  0.26  designates  the  ring  area  where 
the  sail  is  mounted.  These  contour  plots  are  generated  by  using  the  C/  measurements 
obtained  from  both  sensor  sets  A  and  B.  From  figure  5.55,  the  symmetry  in  the  Cf 
distribution  with  respect  to  the  sail  chord  line  at  a  =  0.9  °  can  be  observed.  This 
symmetry  disappears  as  the  angle  of  attack  is  increased.  In  figure  5.62,  at  o  =  15.3  °, 
a  low  velocity  region  can  be  noticed  at  around  x/L  =  0.27  and  (p  =  245  °  which  also 
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matches  with  the  oil  flow  visualization  results  at  that  location  given  by  figure  5.3.  Cf 
on  the  hull  takes  relatively  large  values  at  regions  closer  to  the  sail.  At  higher  angles  of 
attack  starting  from  a  =  21.4  the  increase  in  the  Cj  magnitude  on  the  leeward  side  of 
the  sail  leading  edge  where  the  flow  is  most  accelerated  can  also  be  seen. 


5.2  Steady  Pressure  Measurements 

5.2.1  Data  Reduction 

Pressure  data  reduction  was  straight  forward,  however  some  corrections  to  the  data  had 
to  be  made.  As  described  in  chapter  2,  the  pressure  taps  were  on  the  sensor  plugs  located 
just  behind  each  hot-film  sensor.  Although  the  sensor  plugs  were  levelled  with  the  surface 
of  the  model  as  much  as  possible,  there  were  slight  surface  irregularities  in  the  vicinity 
of  the  pressure  holes  which  effected  the  pressure  distribution.  These  irregularities  were 
minimized  before  the  final  measurements.  However,  there  were  still  some  oscillations 
in  the  data.  In  order  to  get  rid  of  these  oscillations,  the  following  approach  has  been 
developed:  the  local  irregularities  on  the  pressure  ports  effect  the  shape  of  the  local 
streamlines  (make  more  convex  or  concave)  changing  the  local  static  pressure  measured. 
By  using  the  normal-to-wall  momentum  equation,  this  change  in  the  pressure  coefficient 
ACp  can  be  approximated  as: 

(5.5) 

dr  r 

where  U  is  the  velocity  in  the  local  streamwise  direction  and  r  is  the  normal-to-wall 
coordinate.  Assuming  that  the  local  streamline  curvature  due  to  the  surface  irregularity 
will  be  approximately  the  same  at  all  angles  of  attack  and  roll  position  of  the  model; 

Here  A  is  the  unknown  coefficient  which  is  different  for  each  port  position.  With  a  simple 
approach,  Cp  can  be  written  as: 

(5.7) 

^OO 

By  using  the  above  result  in  equation  5.6; 


ACp  =  A{1  -  Cp) 


(5.8) 
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If  we  define  ACp  as  the  correction  to  be  made  to  each  port  location,  then  we  can  write 

C-pc  —  Cpm  +  ACp  (5-9) 

where  Cpm  is  the  measured  Cp  value  and  Cpc  is  the  corrected  Cp  value.  By  using  the 
expression  for  ACp  given  by  equation  5.8,  we  can  obtain  the  final  form  of  the  correction 
equation  as: 

Cpc  =  ^  +  (1  —  A)Cpm  (5.10) 

In  order  to  find  the  A  value  for  each  port  position,  the  Cp  results  of  the  barebody  at  0  ° 
angle  of  attack  obtained  from  a  BANS  code  were  used  (figure  5.69).  These  CFD  results 
were  supplied  by  Dr.  D.  L.  Whitfield,  Mississippi  State  University.  At  each  port  position, 
the  Cpc  was  taken  as  the  Cp  value  of  the  CFD  and  the  Cp  values  of  the  barebody  pressure 
measurements  at  0  °  angle  of  attack  with  slotted  walls  were  taken  as  the  Cpm-  The  A 
values  obtained  as  the  result  of  this  procedure  were  used  for  the  correction  of  the  rest  of 
the  data  that  were  taken  at  10  °  and  20  °  angles  of  attack. 


5.2.2  Results  and  Discussion 

Figures  5.70  and  5.71  show  barebody  pressure  coefficient  Cp  vs.  circumferential  location 
(f)  distributions  for  different  xjL  locations  at  a  =  10  Each  figure  has  10  x/L  locations 
which  gives  a  total  number  of  20  streamwise  measurement  stations  at  a  specific  angle  of 
attack.  The  first  station  is  located  at  x/L  =  0.112  and  the  last  one  at  x/L  =  0.845.  Note 
that  the  pressure  measurement  stations  are  slightly  downstream  of  the  corresponding 
skin-friction  measurement  locations.  However,  the  difference  is  approximately  0.2%  of 
the  overall  model  length  and  can  be  considered  as  negligible.  All  measurement  stations 
given  in  Cp  figures  lie  on  the  long  sensor  row  of  the  model.  Figures  5.72  and  5.73  show 
the  same  kind  of  distribution  for  a  =  20  °. 

In  figures  5.70  to  5.73,  Cp  values  are  presented  for  two  wind  tunnel  wall  configurations. 
Solid  symbols  represent  the  solid  wall  results  and  the  open  symbols  stand  for  the  slotted 
wall  case.  As  can  be  seen  from  these  figures,  Cp  values  of  the  solid  wall  case  are  more 
negative  compared  to  the  ones  obtained  in  the  slotted  wall  case  through  all  </»  locations. 
This  is  an  expected  result,  since  in  the  solid  wall  case  the  fluid  over  the  measurement 
stations  is  accelerated  more  compared  to  the  real  flow  conditions  due  to  the  blockage 
effect  of  the  model  at  angle  of  attack.  This  extra  pressure  gradient  term  makes  the 
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C-p  values  more  negative.  On  the  other  hand,  in  the  slotted  wall  configuration  with  an 
open-air  ratio  of  38%,  the  blockage  effect  can  be  reduced  significantly  and  more  reliable 
pressure  and  skin-friction  values  can  be  obtained.  The  difference  between  the  slotted  and 
solid  wall  results  is  more  significant  for  a  =  20 

At  both  angles  of  attack,  a  favorable  pressure  gradient  can  be  observed  on  the  windward 
side  of  the  model  for  each  sensor.  The  leeward  side  of  model  can  be  thought  as  a  pressure 
recovery  region.  However  due  to  the  crossflow  separation  on  the  leeward  side  of  the  model, 
the  general  Cp  trend  does  not  follow  a  monotonic  increase.  Especially  at  a  =  20  ° ,  a  local 
minimum  in  Cp  distribution  can  be  detected  on  the  leeward  side  of  the  model  at  each 
station  between  x/L  =  0.288  and  x/L  =  0.777.  As  also  observed  by  Wetzel  [2]  in  the 
prolate  spheroid  case,  these  minima  are  the  result  of  the  flow  separation  and  may  coincide 
beneath  the  coherent  vortices  that  form  in  such  an  open  separation.  These  vortices  with 
relatively  high  circulation  induce  again  higher  local  velocities  which  in  turn  lead  to  a  local 
pressure  decrease  immediately  beneath  the  vortical  core.  The  Cp  distributions  of  the  last 
two  stations  x/L  —  0.821  and  x/L  =  0.865  follow  a  different  trend  on  the  leeward  side 
of  the  model  at  both  angles  of  attack  compared  to  the  other  stations.  As  in  the  skin- 
friction  case,  a  relatively  flat  Cp  profile  can  be  observed  on  the  circumferential  locations 
corresponding  to  the  low  speed  separated  flow  region  of  the  stern. 

Figures  5.74  and  5.75  show  Cp  vs.  (j)  distributions  for  different  x/L  locations  at  a  =  10° 
for  the  sail-on-side  case.  The  data  are  compared  with  the  barebody  case  on  the  non-sail 
region  of  the  model.  Solid  symbols  represent  the  Cp  values  of  the  sail-on-side  case  while 
the  open  symbols  are  used  for  the  barebody  case.  Both  sail-on-side  and  barebody  results 
are  obtained  with  the  slotted  wall  configuration.  Figures  5.76  and  5.77  show  the  same 
kind  of  distribution  for  a  =  20  °.  Note  that  the  scales  of  the  Cp  axes  are  different  for 
a  =  10  °  and  a  =  20  °  cases.  On  the  non-sail  region  of  the  model  (between  cf)  =  0° 
and  4>  =  180  °),  the  Cp  distributions  follow  the  same  trend  as  the  one  obtained  for  the 
barebody  case  at  both  angles  of  attack  and  the  values  are  approximately  the  same  within 
the  uncertainties.  On  the  sail  side,  the  mean  pressure  distributions  are  affected  by  the 
presence  of  the  sail  except  the  first  four  stations  on  the  nose  part  of  the  model.  These 
four  stations  show  a  symmetric  pressure  distribution  with  respect  to  4>  —  180  °.  For  the 
other  stations,  this  symmetry  can  not  be  observed.  The  effect  of  the  sail  is  significant  at 
the  stations  located  between  x/L  =  0.210  and  xjL  =  0.308.  The  station  &X  x/L  ~  0.210 
is  just  upstream  of  the  sail  and  the  pressure  increase  due  the  stagnation  process  can  be 


5.2  Steady  Pressure  Measurements 


58 


clearly  seen  between  4>  =  260  °  and  (/>  =  280  The  Cp  value  takes  its  maximum  value 
approximately  at  0  =  280  °  for  both  angles  of  attack.  The  pressure  distributions  in  the 
vicinity  of  the  leading  edge  of  the  sail  are  of  particular  interest,  since  they  are  responsible 
for  the  formation  of  the  horseshoe  vortices  in  this  region  as  described  by  Simpson  [31]. 
During  the  measurements  two  stations  were  under  the  sail:  station  xjL  =  0.268  between 
(j)  =  262  °  and  (j)  =  270  and  station  x/L  =  0.288  between  (p  =  267  °  and  4>  =  269  °.  For 
these  two  stations,  the  Cp  values  at  the  (p  locations  under  the  sail  are  left  blank  in  the 
figures. 

Figures  5.78  through  5.81  show  the  Cp  vs.  xjL  distributions  obtained  with  slotted  walls 
for  different  (p  locations  at  10°  and  20°  angles  of  attack  .  Figures  5.78  and  5.79  give  the 
results  for  the  barebody  case.  For  a  =  10°,  a  smooth  change  of  Cp  with  x/L  at  different 
(p  values  can  be  observed.  At  a  =  20  °,  Cp  values  are  more  scattered  and  the  magnitude 
difference  between  each  (p  position  is  bigger.  Sail-on-side  results  are  given  in  figures  5.80 
and  5.81.  At  both  angles  of  attack,  the  effect  of  the  sail  on  the  pressure  distribution  can 
be  seen  on  the  leeward  side  of  the  model,  especially  at  (/•  =  210  °  and  (p  =  270  °. 

In  order  to  determine  the  circumferential  behavior  of  the  mean  pressure  in  the  separated 
flow  regions,  Cp  distributions  are  compared  with  the  C/  variations  at  the  same  x/L 
locations.  C/  distributions  of  x/L  =  0.501  and  x/L  =  0.774  at  a  =  11.3°  and  a  =  21.4° 
are  compared  with  the  Cp  variations  of  x/L  =  0.503  and  x/L  =  0.777  at  a  =  10  ° 
and  a  =  20  °.  Figures  5.82  to  5.85  show  the  comparison  of  the  barebody  case  at  the 
conditions  given  above,  while  the  figures  5.86  to  5.89  show  the  comparison  of  the  sail- 
on-side  case.  A  common  characteristic  of  the  pressure  distributions  in  the  vicinity  of 
the  separation  locations  can  be  observed  from  the  barebody  figures:  the  mean  pressure 
values  are  approximately  constant  over  the  separated  fluid  regions.  This  zero  pressure 
gradient  region  is  obvious  at  x/L  =  0.777  for  a  =  20  °  given  by  the  figure  5.85.  The  flat 
pressure  profile  on  the  leeward  side  of  the  model  covers  the  region  between  (p  =  111  ° 
and  (p  =  153  °  and  includes  both  the  primary  and  the  secondary  separation  locations. 
The  same  Cp  trend  can  be  observed  on  the  leeward  side  of  the  non-sail  region  in  the 
sail-on-side  case.  As  discussed  before,  at  relatively  low  angles  of  attack,  two  separation 
locations  can  be  observed  on  the  sail  region  of  model.  In  figure  5.86,  these  two  separation 
locations  can  be  seen  on  the  leeward  side  of  the  sail  region  at  x/L  =  0.501  for  a  =  11.3°. 
Over  the  first  separation  location  {(p  =  216  °),  the  pressure  gradient  is  not  zero  whereas 
for  the  second  separation  location  {(p  =  260°)  a  constant  pressure  region  exists.  At  higher 
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angles  of  attack,  only  the  second  separation  location  remains  and  the  first  one  vanishes. 
In  figure  5.89,  at  x/L  =  0.777  for  a  =  20°,  again  a  flat  pressure  profile  over  the  separation 
location  on  the  leeward  side  of  the  sail  region  can  be  detected.  As  discussed  in  Wetzel 
et  al.  [32],  although  the  pressure  data  are  used  often  to  indicate  the  existence  of  massive 
separation,  it  is  not  a  good  indicator  of  the  separation  location.  This  is  due  to  the  fact 
that  the  crossflow  separation  is  the  result  of  a  local  flow  phenomena,  but  pressure  at  a 
given  point  in  space  is  strongly  influenced  by  the  entire  flow  field.  The  results  of  this 
study  also  support  the  aforementioned  principle.  One  can  think  of  using  the  flat  pressure 
profile  in  regions  of  separation  to  locate  the  separation  locations,  however  it  is  difficult 
to  determine  the  exact  point  where  this  flat  pressure  distribution  begins. 
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Figure  5.1:  Oil  flow  visualization  showing  the  cross  flow  separation  topology  on  the  constant 
diameter  region  of  the  model  for  the  barebody  case  at  a  =  20  Re  =  4.5  x  10®.  Flow  is 
from  left  to  right. 


Figure  5.2:  Oil  flow  visualization  showing  the  low  speed  fluid  region  on  the  stern  of  the  model 
for  the  barebody  case  at  a  =  15  Re  =  4.5  x  10®.  Flow  is  from  left  to  right. 
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Figure  5.3:  Oil  flow  pattern  showing  the  separation  in  the  vicinity  of  the  sail  at  a  =  15 
Rr  =  4.5  X  10®.  Flow  is  from  right  to  left. 


Figure  5.4:  Oil  flow  visualization  of  the  sail  side  showing  the  low  speed  fluid  region  on  the 
stern  of  the  model  for  the  sail-on-side  case  at  a  =  10  Re.  =  4.5  x  10®.  Flow  is  from  right 
to  left. 
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Figure  5.8:  Cf  vs.  </>  for  all  x/L  locations  at  a  =  7.2  °  for  steady  barebody  and  sail-on-side 
cases.  Sail  side  on  the  right  of  the  figure  starting  from  cj)  =  180  . 
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Figure  5.9:  C/  vs.  4>  for  all  x/L  locations  at  a  =  9.3°  for  steady  barebody  and  sail-on-side 
cases.  Sail  side  on  the  right  of  the  figure  starting  from  (j)  =  180  °. 
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Figure  5.12:  Cf  vs.  (p  for  all  x/L  locations  at  a  =  15.3°  for  steady  barebody  and  sail-on-side 
cases.  Sail  side  on  the  right  of  the  figure  starting  from  (p  =  180  °. 
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Figure  5.13;  C/  vs.  <j)  for  all  x/L  locations  at  a  =  17.4°  for  steady  barebody  and  sail-on-side 
cases.  Sail  side  on  the  right  of  the  figure  starting  from  (f>  =  180  °. 
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Figure  5.18:  Cf  vs.  <p  for  all  xjL  locations  at  a  =  27.6°  for  steady  barebody  and  sail-on-side 
cases.  Sail  side  on  the  right  of  the  figure  starting  from  (f)  =  180  °. 
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Figure  5.19:  Cf  vs.  0  for  all  a  at  x/L  =  0.110  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  (j)  =  180  °. 
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Figure  5.20;  Cj  vs.  0  for  all  a  at  x/L  =  0.131  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  ^  =  180  °. 
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Figure  5.21;  Cf  vs.  0  for  all  a  at  x/L  =  0.170  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  (j)  =  180  °. 
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Figure  5.22:  Cf  vs.  0  for  all  a  at  x/L  =  0.266  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  0  =  180  °. 
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Figure  5.24:  C/  vs.  <p  for  all  a  at  x/L  =  0.345  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  (f)  =  180  °. 
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Figure  5.25:  C/  vs.  <p  for  all  a  at  x/L  =  0.434  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  <j)  =  180  °. 
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Figure  5.26:  Cf  vs.  </>  for  all  a  x/L  —  0.501  for  steady  barebody  and  sail-on-side  cases 
Sail  side  on  the  right  of  the  figure  starting  from  (f)  =  180  °. 
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Figure  5.29:  Cf  vs.  (j)  for  all  a  at  x/L  =  0.706  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  (/>  =  180  °. 
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Figure  5.30:  C/  vs.  (j)  for  all  a  at  x/L  =  0.774  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  0  =  180  °. 
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Figure  5.32;  C/  vs.  (f>  for  all  a  at  xjL  =  0.863  for  steady  barebody  and  sail-on-side  cases. 
Sail  side  on  the  right  of  the  figure  starting  from  (f)  =  180  °. 
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Figure  5.34:  Steady  primary  separation  locations  vs.  a  for  the  barebody  case  (continued) 
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Figure  5.35:  Steady  secondary  separation  locations  vs.  a  for  the  barebody  case. 


Figure  5.36:  Steady  secondary  separation  locations  vs.  a  for  the  barebody  case  (continued). 
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Figure  5.37;  Comparison  of  the  oil-flow  primary  separation  locations  with  the  hot-film  primary 
separation  locations  for  the  barebody  case.  For  the  hot-film  measurements  a  =  15.3  °  and 
for  the  oil-flow  results  a  =  15  °. 


Figure  5.38:  Comparison  of  the  oil-flow  primary  separation  locations  with  the  hot-film  primary 
separation  locations  for  the  barebody  case.  For  the  hot-film  measurements  a  =  19.4  °  and 
for  the  oil-flow  results  a  =  20  °. 
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Figure  5.42;  Solid  vs.  slotted  wall  comparison  for  the  barebody  Cf  vs  cf)  distribution  at 
xjL  =  0.570.  a  =  20  °  for  the  solid  wall  case,  and  a  —  19.4  °  for  the  slotted  wall  case. 
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Figure  5.45:  Steady  primary  separation  locations  vs.  a  for  the  sail-on-side  (region  without 
sail)  case. 
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Figure  5.46;  Steady  primary  separation  locations  vs.  a  for  the  sail-on-side  (region  without 
sail)  case  (continued). 
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Figure  5.49:  Comparison  for  steady  primary  separation  locations  vs.  a  for  barebody  and  the 
sail-on-side  (region  without  sail)  case  at  x/L  =  0.501,  0.638  and  0.774. 
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Figure  5.50:  Comparison  for  steady  secondary  separation  locations  vs.  a  for  barebody  and 
the  sail-on-side  (region  without  sail)  case  at  x/L  =  0.501,  0.638  and  0.774. 
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Figure  5.52:  Steady  1"*  separation 
(continued). 
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Figure  5.54:  Steady  2^  separation  locations  vs.  a  for  the  sail-on-side  (region  with  sail)  case 
(continued) 


Figure  5.55:  Steady  Cf  contours  in  the  vicinity  of  sail  at  a  =  0.9 


Figure  5.56:  Steady  Cr  contours  in  the  vicinity  of  sail  at  a  =  3.1 
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Figure  5.65:  Steady  Cj  contours  in  the  vicinity  of  sail  at  a  =  21.4  °. 
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Figure  5.66:  Steady  Cf  contours  in  the  vicinity  of  sail  at  a  =  23.4  °. 
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Figure  5.67:  Steady  Cj  contours  in  the  vicinity  of  sail  at  a  =  25.5  °. 
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Figure  5.68:  Steady  (7/  contours  in  the  vicinity  of  sail  at  a  =  27.6  °. 
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Figure  5.70:  Cp  vs.  (j)  for  the  stations  between  x/L  —  0.112  and  xjL  —  0.327  at  a  —  10.0 ' 
for  the  barebody  case.  Solid  and  slotted  wind  tunnel  wall  configurations. 
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Figure  5.71:  Cp  vs.  (/>  for  the  stations  between  x/L  =  0.348  and  x/L  =  0.865  at  a  =  10. 
for  the  barebody  case.  Solid  and  slotted  wind  tunnel  wall  configurations. 
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Figure  5.72:  Cp  vs.  ^  for  the  stations  between  xjL  =  0.112  and  x/L  =  0.327  at  a  =  20.0" 
for  the  barebody  case.  Solid  and  slotted  wind  tunnel  wall  configurations. 
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Figure  5.73;  Cp  vs.  (f)  for  the  stations  between  x/L  =  0.348  and  xjL  =  0.865  at  a  =  20. 
for  the  barebody  case.  Solid  and  slotted  wind  tunnel  wall  configurations. 
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Figure  5.74:  Cp  vs.  0  for  the  stations  between  xjL  —  0.112  and  xjL  =  0.327  at  a  =  10.0° 
for  the  barebody  case  and  sail-on-side  cases.  Sail  side  on  the  right  of  the  figure  starting  from 
(f)  =  180  °.  Slotted  wind  tunnel  wall  configuration. 
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Figure  5.75:  Cp  vs.  <j)  f< 
for  the  barebody  case  ai 
4>  =  180°.  Slotted  wine 
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Figure  5.77:  Cp  vs.  (p  for  the  stations  between  x/L  =  0.348  and  x/L  =  0.865  at  a  =  20.0° 
for  the  barebody  case  and  sail-on-side  cases.  Sail  side  on  the  right  of  the  figure  starting  from 
</)  =  180  °.  Slotted  wind  tunnel  wall  configuration. 
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Figure  5.78:  Barebody  Cp  vs.  x/L  distributions  for  different  cj)  locations  at  a  —  10  °. 


Figure  5.79:  Barebody  Cp  vs.  x/L  distributions  for  different  (j)  locations  at  a  —  20  °. 
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Figure  5.80:  Sail-on-side  C-p  vs.  x/L  distributions  for  different  0  locations  at  a  =  10  °. 


Q- 

o 


0.30 

0.20 

0.10 

0.00 

-0.10 

-0.20 

-0.30 

-0.40 

-0.50 

-0.60 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 

X/L 


0.6 


0.7 


0.8 


0.9 


a  =  20^ 

□  phi  =  0  deg. 

A  phi  =  180  deg. 

•  phi  =  210  deg. 

0  phi  =  240  deg. 
+  phi  =  270  deg 

- Bn - 

□nnrn'-'°  B  □  n  B 

> 

«0 

• 

□ 

— - - ^ 

• 

A  f  6  ♦  o 

AA-  aA  _i_  +  ▼  1 

0 

- JU - 

-I-++  5  A  ^ 

^  .1-  •  ^ _ A - 

D  5  □ 

0  D 

+ 

• 

0 

> 

1 

> 

? 

- A: — - 

_ _ _ 

4- 

,  - 1 - -1 - 1 . . -n - “-I - -  I  '  . '  “ 

1.0 


Figure  5.81;  Sail-on-side  Cp  vs.  xjL  distributions  for  different  0  locations  at  a  —  20°. 
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Figure  5.82:  Steady  Cp  and  Cf  vs.  (j)  for  the  barebody  case.  For  Cp  measurements  a  =  10.0° 
and  x/L  =  0.503.  For  C/  measurements  a  =  11.3°  and  xlL  =  0.501. 


Figure  5.83;  Steady  Cp  and  Cj  vs.  (f)  for  the  barebody  case.  For  Cp  measurements  a  =  20.0 
and  xjL  =  0.503.  For  C/  measurements  a  —  21.4°  and  x/L  =  0.501. 
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Figure  5.84:  Steady  Cp  and  C/  vs.  </>  for  the  barebody  case.  For  Cp  measurements  a  —  10.0° 
and  xjL  =  0.777.  For  Cj  measurements  a  -  11.3°  and  x/L  =  0.774. 


Figure  5.85:  Steady  Cp  and  Cj  vs.  0  for  the  barebody  case.  For  Cp  measurements  a  -  20.0° 
and  x/L  =  0.777.  For  Cf  measurements  a  =  21.4°  and  x/L  =  0.774. 


Chapter  6 


Results  and  Discussion  for  Unsteady 
Measurements 


6.1  Data  Acquisition  and  Reduction 

The  unsteady  volta,ge  va.lues  from  the  constant  temperature  anemometers,  each 

connected  to  a  hot-film  sensor,  were  acquired  by  using  a  custom  written  Labview  software 
program.  This  program  also  provided  the  control  and  the  adjustment  of  the  tunnel  data 
acquisition  hardware  for  the  unsteady  hot-film  measurements.  The  number  of  samples 
taken  in  time  for  each  unsteady  maneuver  was  1500.  The  unsteady  voltage  signals  were 
acquired  with  a  sample  rate  of  1000  Hz.  Note  that  sampling  frequency  is  well  above  the 
frequency  response  of  the  constant  temperature  anemometers  (200  Hz).  Total  sampling 
time  was  1.5  seconds.  The  actual  maneuver  covered  approximately  0.33  seconds  of  this 
time  period.  The  rest  of  the  sampling  time  was  spent  for  taking  the  pre-maneuver  and 
the  post-maneuver  data.  The  pitch-up  maneuver  was  executed  10  times  at  each  roll  angle 
position  in  order  to  reduce  the  uncertainty  of  unsteady  data. 

As  used  in  the  steady  data  reduction  part,  the  measured  unsteady  values  are  represented 
as  arrays  of  indexed  variables  in  order  to  simplify  the  description  of  the  data  reduction 
procedure.  In  the  unsteady  data  reduction  process,  the  unsteady  voltage  value  acquired 
from  each  anemometer  can  be  written  as 


E{t')  =  E{tj,Sk,ri,(j)m) 


(6.1) 
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j  =  l,...,1500  /c  =  l,...,15  Z  =  l,...,10  and  m  -  1,  ...,nro/Z  (6.2) 

where  tj  corresponds  to  a  point  in  time;  Sk  to  the  sensor  number  or  the  location;  n  to 
the  repetition  number  of  the  pitch-up  maneuver;  and  to  the  DyPPiR  roll  angle  of  the 
model.  Total  number  of  the  roll  angles  covered  during  the  experiments,  nroU,  changes 
with  respect  to  the  model  configuration.  For  the  barebody  case,  56  roll  angle  positions 
were  covered.  For  the  sail-on-side  case,  total  number  of  roll  angles  was  117.  However, 
data  of  the  non-sail  region  and  the  data  from  sail  region  were  reduced  separately.  For 
the  non-sail  region,  total  number  of  71  roll  angles  were  reduced,  while  for  the  sail  side 
data  reduction  was  made  with  76  roll  angles.  The  last  30  angles  of  the  non-sail  region 
and  the  first  30  angles  of  the  sail  side  were  the  same.  This  was  done  in  order  to  keep  the 
continuity  of  the  skin-friction  distributions  between  the  two  regions  during  the  smoothing 
process,  which  will  be  described  later  in  this  section. 

The  steps  of  the  unsteady  skin-friction  data  reduction  can  be  outlined  as  follows: 


1.  At  each  roll  angle  position,  the  free-stream  temperature  T{tj,ri,(f)m),  the  dynamic 
pressure  (pm),  and  the  atmospheric  pressure  p{tj,ri,  <pm)  are  calculated  from 

the  corresponding  voltage  values  by  the  use  of  the  calibration. 


2.  At  each  roll  angle  position,  the  plunge  position  z{tj,ri,(pm)  and  the  instantaneous 
angle  of  attack  a{tj,ri,(pm)  of  the  model  are  determined  through  the  use  of  the 
DyPPiR  plunge  and  pitch  calibrations. 


3.  Each  time  point  is  non-dimensionalized  by  the  reference  time,  tref  =  LjU^  where 
L  is  the  model  length  and  Uoc  is  the  free-stream  velocity. 

4.  Instantaneous  skin  friction  values  Cf{tj,Sk,ri,(pm)  for  each  roll  angle,  sensor  and 
the  repetition  are  calculated  by 


Sk,  r/,  ^m) 


At) (t j ,  Sk,f'l,  (pm) 

qi^j,  I'l,  ^m) 


(6.3) 


where  is  the  unsteady  wall  shear  stress  value  obtained  by  the  use  of  King’s  law 
given  by  the  equation  4.10 


At)(^j)  ^k,  D)  ^m) 


Sk,  ri,<pm)  -  A{Sk,  (pm) 
B{Sk,^m) 


3 


(6.4) 
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In  equation  6.4,  A{sk,4>m)  the  calibration  coefficients  obtained 

for  each  sensor  by  the  procedure  described  in  chapter  IV.  Note  that  at  a  particular 
roll  angle  position  of  the  model,  the  calibration  coefficients  used  in  the  unsteady 
measurements  are  the  same  as  the  ones  used  for  the  steady  measurements,  since  the 
maneuvers  are  performed  just  after  the  last  steady  angle  of  attack  and  the  change 
in  the  free-stream  temperature  between  the  steady  and  the  unsteady  measurements 
at  a  specific  roll  angle  is  negligible. 

5.  To  reduce  the  noise  in  the  unsteady  measurements,  the  skin-friction  values  of  each 
sensor,  the  plunge  position  and  instantaneous  angle  of  attack  from  10  repetitions 
of  the  pitch-up  maneuver  performed  at  a  specific  roll  angle  position  of  the  model 
are  ensemble  averaged: 


Cfa„,{tj,Sk,(l>m}  -  jQ 

(6.5) 

ZavgV'jiym)  —  jq 

(6.6) 

(6.7) 

6.  Standard  deviations  a{tj,Sk,<l>m)  of  the  ensemble  averaged  skin-friction  values  are 
calculated.  The  results  from  the  unsteady  experiments  showed  that  the  standard 
deviations  of  the  ensemble  averaged  voltage  values  were  two  orders  of  magnitude 
smaller  than  the  average  values  in  most  cases,  therefore  the  unsteady  data  was 
highly  repeatable  and  10  repetitions  for  the  pitchup  maneuver  at  each  roll  angle 
position  were  sufficient  to  reduce  the  noise. 

7.  In  order  to  determine  the  bad  skin-friction  values  and  exclude  them  from  ensem¬ 
ble  averaging,  a  statistical  outlier  detection  method  is  used.  In  order  to  deter¬ 
mine  the  outliers,  a  student's  t  -  distribution  is  assumed  for  the  10  repetitions 
of  the  skin-friction  values  Cf{tj,Sk,ri,(f>rn)  at  each  roll  angle  position.  Then  the 
difference  d{tj,Sk,<f>Tn)  from  the  mean  value  Cf^„g{tj,  Sk,(l>m)  is  calculated  for  each 

d{tj,Sk,<pm)  —  \J [Cfi'tjt  SkiTl,(f)m)  ~  */*"»)]  (®'^) 

and  the  condition 

d>t  X  a{tj,Sk,  4>m)  (6-®) 
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is  checked.  Here  t  value  is  chosen  for  9  degrees  of  freedom  with  90%  confidence 
interval  and  equal  to  1.833  (Holman  [33]).  The  Cf{tj,Sk,riAm)  values  that  satisfy 
the  above  condition  are  marked  as  outliers  and  their  values  are  set  to  zero.  The 
new  are  calculated  by: 


TJL  Cf{tj,Sk,ri,<t>m) 

10  '^outlier 


(6.10) 


where  TioutUer  is  the  total  number  of  outliers  detected. 


8.  In  order  to  reject  the  high  frequency  noise  in  Sfc,</>m),  a  third-order  digital 

Butterworth  filter  is  used.  As  for  the  cutoff  frequency,  20  Hz.  is  used.  Wetzel  [2] 
used  the  same  value  in  the  unsteady  data  reduction  of  the  prolate  spheroid  work. 
In  this  study,  different  cutoff  frequency  values  were  tried.  It  has  been  seen  that 
this  value  is  high  enough  to  preserve  the  frequency  content  and  the  magnitude  of 
the  unsteady  data  while  filtering  the  high  frequency  noise. 


9.  Unsteady  separation  locations  are  obtained  and  the  data  are  arranged  for  different 
output  formats. 


6.1.1  Determination  of  the  unsteady  separation  locations 

For  the  determination  of  the  unsteady  primary  separation  locations,  Cf  =  C/„„^  {tj  ,Sk,<t>m) 
vs.  (^  =  (^  =  1)  ■■■inroll)  distributions  for  each  sensor  have  to  be  analyzed.  Since  the 

actual  maneuver  starts  at  t'  =  3.0  and  ends  at  approximately  f  =  10.3,  there  is  no  need 
to  consider  all  the  time  steps  in  the  maneuver.  Between  these  two  time  values,  a  total 
number  of  40  equally  spaced  time  points  are  used.  This  corresponds  to  a  time  increment 
of  At  =  0.01  seconds  between  the  start  of  the  pitch-up  and  the  end  of  the  maneuver. 
The  total  number  of  time  points  to  be  used  in  the  separation  location  determination 
is  obtained  as  a  result  of  many  trials  and  it  has  been  seen  that  further  increase  of  the 
total  number  of  points  does  not  change  the  unsteady  separation  pattern  obtained  with 
40  points. 

An  automated  scheme  is  needed  to  find  the  unsteady  separation  locations  at  each  time 
step  and  the  sensor  location,  since  it  is  impossible  to  determine  the  separation  locations 
graphically  (as  done  in  the  steady  case)  due  to  the  large  number  of  data  sets  to  be  exam¬ 
ined.  For  this  purpose,  three  methods  were  tried:  in  the  first  method,  the  minimum  in 
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Cf  from  each  unsteady  circumferential  distribution  is  determined  and  the  corresponding 
(f)  is  taken  to  be  the  unsteady  separation  location.  However,  this  method  is  not  robust 
and  accurate  enough  in  the  presence  of  outliers  or  wiggles  in  the  separation  region  and 
gives  misleading  results.  In  the  second  method,  after  the  minimum  in  (7/  is  determined, 
a  quadratic  curve  is  fit  to  the  concave  separation  region  by  the  least  squares  method.  In 
the  fit  process  11  (C/,0)  points  are  used  with  the  center  value  being  the  minimum  of  C/ 
obtained  from  the  first  method.  This  helps  smoothing  the  data  near  the  separation  region 
and  improves  the  separation  location  determination.  However,  as  observed  in  most  of  the 
actual  cases,  the  Cf  distribution  in  the  separation  region  is  not  symmetric  with  respect  to 
the  minimum  Cf  value,  therefore  using  a  quadratic  fit  shifts  the  minimum  from  the  true 
value.  This  offset  is  not  constant  for  each  distribution,  so  a  bias  correction  to  the  data 
cannot  be  applied.  The  third  method  removed  the  problems  encountered  in  the  previous 
methods.  In  this  approach,  the  loess  method  is  used  as  described  in  the  steady  data 
reduction  of  the  sail-on-side  case.  This  method  is  considered  in  order  to  smooth  the  Cf 
vs.  0  distributions  at  each  time  step  without  changing  the  circumferential  skin-friction 
pattern  and  to  exclude  the  deviant  points  that  distort  the  smoothed  data.  By  using  this 
technique,  the  asymmetry  in  the  vicinity  of  the  separation  location  is  preserved,  while 
the  noise  from  the  data  is  removed.  The  magnitude  changes  in  the  skin-friction  values 
due  to  smoothing  is  found  to  be  within  the  uncertainties.  As  a  result,  this  last  method  is 
used  to  find  the  unsteady  separation  locations  and  to  obtain  the  Cf  vs.  0  distributions 
at  each  time  step. 


6.2  Unsteady  Flow  Topology 

One  of  the  methods  to  investigate  the  difference  between  the  steady  and  the  unsteady 
flow  over  the  model  is  to  compare  the  steady  (7/  vs.  0  distribution  obtained  at  a  certain 
pitch  angle  a  with  the  unsteady  one  acquired  at  the  corresponding  instantaneous  pitch 
angle  a{t')  for  a  specific  xjL  location.  Figures  6.1  through  6.39  show  this  comparison 
for  all  the  xjL  locations  of  the  sensor  set  B.  Figures  from  6.1  to  6.13  give  the  barebody 
results.  Figures  6.14  to  6.26  show  the  data  of  the  non-sail  region  of  the  sail-on-side  case, 
and  the  figures  between  6.27  and  6.39  give  the  comparison  of  the  region  with  the  sail. 
Note  that,  for  each  figure,  the  steady  angle  of  attack  and  the  instantaneous  angle  of 
attack  are  the  same.  The  non-dimensional  time  t'  corresponding  to  the  instantaneous 
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angle  of  attack  is  also  shown  in  each  figure.  All  the  unsteady  experiments  were  done  with 
the  slotted  walls. 

Unsteady  Cf  distribution  of  the  barebody  case  and  the  non-sail  region  of  the  sail-on- 
side  case  show  similar  trends  compared  to  the  steady  distributions  of  each.  As  seen 
from  figures  6.1  and  6.14,  at  a  =  3.1  the  skin-friction  magnitudes  of  the  steady  and 
the  unsteady  cases  are  approximately  the  same  and  the  distributions  of  both  are  flat 
profiles.  As  the  steady  angle  of  attack  is  increased,  the  steady  data  start  to  deviate  from 
this  flat  profile,  and  minima  in  Cf  distributions  indicating  primary  separation  can  be 
observed.  However,  the  unsteady  skin-friction  pattern  remains  approximately  flat  until 
a(t')  =  13.2  °.  At  this  instantaneous  angle  of  attack,  the  unsteady  primary  separation 
starts  to  develop  at  the  downstream  stations  starting  from  x/L  =  0.706.  Figures  6.6 
and  6.19  show  the  Cf  distributions  at  this  angle  of  attack  for  the  barebody  and  the 
sail-on-side  case  respectively.  At  higher  angles  of  attack,  the  difference  between  the 
steady  and  the  unsteady  flow  separation  pattern  is  more  obvious.  At  a  =  21.4  °,  in 
figures  6.10  and  6.23,  steady  primary  separation  starts  from  x/L  =  0.266  and  shifts 
towards  the  windward  side  at  the  downstream  stations.  The  unsteady  separation  can  be 
first  located  at  x/L  =  0.345  and  also  moves  towards  the  windward  side  as  x/L  increases. 
Compared  to  the  steady  primary  separation  locations  at  the  same  x/L  stations,  the 
unsteady  separation  locations  are  found  to  be  more  leeward.  In  particular,  ai  x/L  ■= 
0.501,  the  primary  separation  location  for  the  steady  data  is  0  =  122°,  while  the  unsteady 
primary  separation  occurs  aX  4>  =  136  °  for  the  barebody  case.  From  the  same  figures, 
another  important  observation  can  be  made  about  the  secondary  flow  separation.  In  the 
steady  case,  secondary  separation  starts  to  develop  at  x/L  =  0.434,  and  can  be  seen 
as  the  second  minima  of  the  Cf  for  the  rest  of  the  stations  downstream.  On  the  other 
hand,  no  secondary  separation  formation  can  be  detected  at  any  of  the  x/L  stations 
for  the  unsteady  case.  A  weak  unsteady  secondary  separation  can  be  seen  at  the  last 
instantaneous  angle  of  attack  a(t')  =  27.6°,  at  the  last  four  x/L  stations.  As  previously 
described  in  the  steady  results,  at  the  first  three  stations  on  the  nose  region,  a  weak 
separation/reattachment  pattern  can  be  seen  on  the  leeward  side  of  the  model  at  high 
steady  angles  of  attack.  However,  this  flow  behavior  cannot  be  observed  for  the  unsteady 
case  at  the  corresponding  instantaneous  angles  of  attack.  In  general,  the  magnitudes  of 
the  unsteady  Cf  are  lower  than  that  of  the  steady  ones  at  the  same  angle  of  attack  and 
the  x/L  location.  This  difference  can  be  observed  starting  from  a  =  5.1  °  and  become 
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more  obvious  at  higher  angles  of  attack,  especially  on  the  leeward  side  of  the  model. 

For  the  sail  region,  the  unsteady  data  start  to  deviate  from  the  steady  Cf  values  even 
at  the  first  instantaneous  angle  of  attack  oi{t')  =  3.1  In  figure  6.27,  the  disturbances 
created  by  the  sail  can  be  seen  on  the  Cf  distribution  at  the  stations  just  downstream 
of  the  sail  for  the  steady  case,  while  the  unsteady  data  do  not  show  the  same  trend 
but  follow  a  more  fiat  pattern.  (Note  that  the  unsteady  Cf  values  at  x/L  =  0.266  are 
not  shown  for  the  <f)  values  after  256  in  figures  6.27  through  6.39,  since  the  data  at 
these  locations  were  too  scattered  and  could  not  be  distinguished  from  the  noise.)  No 
clear  unsteady  separation  pattern  like  the  first  or  the  second  minima  seen  in  the  steady 
case  can  be  observed  up  to  a{t')  =  13.2  °.  After  this  angle  of  attack,  minima  in  the 
Cf  indicating  the  unsteady  separation  can  be  located.  As  also  seen  in  the  steady  case, 
the  unsteady  separation  structure  of  the  sail  side  is  different  from  the  unsteady  crossflow 
separation  topology  observed  on  the  barebody  and  the  non-sail  region  of  the  sail-on-side 
case.  Again,  no  clear  categorization  of  the  separation  locations  could  be  made.  From 
figure  6.38,  for  a{t')  =  25.5°,  the  second  minima  in  the  unsteady  Cf  distributions  at  the 
sensor  stations  downstream  of  xjL  =  0.638  can  be  seen  to  appear  more  leeward  compared 
to  the  second  minima  of  the  steady  case. 

The  results  above  show  that  in  unsteady  flows  the  separation  topology  is  different  from 
the  one  in  an  equivalent  steady  configuration.  For  the  barebody  and  the  non-sail  region 
of  the  sail-on-side  case,  this  difference  originates  from  the  fact  that  the  unsteady  separa¬ 
tion  location  lags  the  steady  separation.  This  can  be  clearly  seen  from  the  figures  6.40 
through  6.59  where  the  quasi-steady  and  the  unsteady  primary  separation  locations  are 
plotted  against  the  non-dimensional  time  t'  for  different  x/ L  stations.  The  quasi-steady 
data  distribution  are  formed  by  first  taking  the  steady  separation  locations  and  the  cor¬ 
responding  steady  angles  of  attack  at  a  certain  x/L  location.  By  using  the  DyPPiR  pitch 
feedback  data,  the  instantaneous  angles  of  attack  that  match  with  the  steady  angles  are 
selected.  Then  the  non-dimensional  time  t'  for  each  steady  separation  location  can  be 
obtained  and  quasi-steady  data  as  a  function  of  t'  can  be  generated.  Open  symbols  stand 
for  the  quasi-steady  data  and  the  filled  symbols  represent  the  unsteady  data.  For  a  bet¬ 
ter  understanding  of  the  mathematical  difference  between  the  steady,  quasi-steady,  and 
the  unsteady  separation  locations,  the  functional  relations  that  are  defined  in  chapter  1 
should  be  referred. 

Figures  6.40  to  6.49  show  the  separation  locations  for  the  barebody  case,  and  the  figures 
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from  6.50  to  6.59  give  the  primary  separation  locations  for  the  non-sail  region  of  the 
sail-on-side  case.  In  these  figures,  the  horizontal  distance  between  the  unsteady  and 
the  steady  separation  location  is  an  indication  of  the  presence  of  a  time  lag.  Since  the 
maneuvers  were  performed  with  a  constant  pitch  rate,  the  instantaneous  pitch  angle  is  a 
linear  function  of  t'.  By  considering  this  fact  in  examining  figures  6.1  to  6.26,  it  can  also 
be  thought  that  at  an  instantaneous  pitch  angle,  for  a  specific  x/L  location,  the  primary 
unsteady  separation  starts  more  leeward  compared  with  the  steady  case.  Therefore  the 
total  separated  region  occupies  a  smaller  area,  both  on  the  model  surface  and  in  the  flow 
field. 

To  understand  why  a  lag  exists,  as  described  by  Wetzel  and  Simpson  [3],  one  should 
examine  the  path  of  a  fluid  particle  during  the  maneuver.  During  the  pitch-up  maneuver, 
a  fluid  particle  originally  starts  its  trajectory  on  the  windward  side  of  the  model  at  a 
lower  angle  of  attack.  Therefore,  the  particle  experiences  a  less  severe  adverse  pressure 
gradient  along  its  path,  and  can  travel  farther  around  the  leeside  before  separating  when 
compared  with  a  similar  particle  in  a  steady  flow  field.  This  history  effect,  the  true  total 
path  the  particle  has  traversed,  differentiates  the  unsteady  flow  field  from  the  steady  one. 

Figures  6.60  through  6.66  show  the  quasi-steady  and  the  unsteady  separation  locations 
as  a  function  of  t'  at  seven  xjL  stations  for  the  sail  region  of  the  sail-on-side  case.  The 
separation  locations  correspond  to  the  second  minima  (at  high  a{t')  these  become  the 
only  minima  as  in  the  steady  case)  of  the  Cf  distributions  measured  from  0  =  180  . 
Figures  6.60  and  6.61  show  that  the  unsteady  separation  locations  are  approximately  the 
same  as  the  steady  ones  within  the  uncertainties  at  x/L  =  0.434  and  xjL  =  0.501.  At 
these  stations,  the  unsteady  separation  locations  are  detected  for  the  time  points  starting 
from  t'  ~  8.0  and  the  </>  values  of  the  separation  locations  are  constant.  At  the  other 
stations,  x/L  =  0.570,  0.638,  0.706,  0.774,  and  0.819,  the  unsteady  separation  locations 
are  more  leeward  compared  to  the  quasi-steady  ones  at  a  given  t'  during  the  maneuver 
and  they  approach  the  steady  values  at  the  end  of  the  maneuver.  The  difference  between 
the  unsteady  and  the  quasi-steady  data  gets  bigger  as  x/L  increases. 

Maybe  the  most  important  observation  that  can  be  made  from  these  figures  is  about 
the  unsteady  separation  pattern  in  the  sail  region  as  a  function  of  t' .  Except  for  x/L  = 
0.434  and  0.501,  the  unsteady  separation  data  do  not  reach  any  of  the  steady  separation 
locations  at  any  instant  of  time  during  the  maneuver.  In  other  words,  unsteady  separation 
pattern  do  not  follow  the  quasi-steady  data  with  a  time  lag.  In  fact,  a  lag  definition  as 
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described  for  the  unsteady  crossflow  separation  may  not  be  appropriate  for  the  sail  side 
separation  topology  based  on  the  results  obtained  in  this  study.  The  results  also  indicate 
the  complex  nature  of  the  unsteady  flow  separation  on  the  sail  side. 


6.3  Time-Lag  Models 

The  main  interest  is  to  be  able  to  model  or  approximate  the  time  lags  associated  with 
the  unsteady  flow  fields.  If  one  can  determine  the  real  physics  behind  the  time-lag 
formation  and  approximate  the  lags  with  reasonable  accuracy,  then  by  using  the  quasi¬ 
steady  data  that  are  constructed  through  the  use  of  the  steady  flow  field  information, 
the  unsteady  flow  features  such  as  unsteady  skin-friction,  force  and  moment  coefficients 
can  be  estimated.  In  the  next  two  sections,  algebraic  time-lag  modelling  and  a  first  order 
differential  time-lag  model  will  be  described  with  their  application  to  the  current  study. 

6.3.1  Algebraic  Time-lag  Models 

Algebraic  time  lag  modeling  was  a  commonly  used  approach  in  the  previous  studies. 
Most  of  the  algebraic  time  lag  models  are  based  on  finding  an  effective  angle  of  attack 
ct^ff  by  using  the  descriptions  of  the  ffow  kinematics.  If  a  model  rotates  about  some  point 
Xcg  at  a  constant  pitch  rate  d,  a  relative  velocity  normal  to  the  model  axis  d(a;  —  Xcg)  will 
exist  at  other  xjL  locations  in  addition  to  the  velocity  of  the  wind  relative  to  the  point 
of  rotation.  By  using  this  fact,  Montividas  et  al.  [9]  approximated  an  effective  angle  of 
attack  as: 

Oiejf  =  a  —  ^OLeff  (611) 

where 

(6.12) 

U  oo 

However,  this  approximation  did  not  describe  the  unsteady  effects  in  their  flow  field 
associated  with  the  onset  of  asymmetric  vortex  shedding  at  high  angles  of  attack  on 
ogive  cylinders. 

Ericsson  [34]  also  studied  the  same  problem  and  extended  the  effective  angle  of  attack 
idea  by  including  the  convective  lag  effects  in  the  flow  field.  He  used  the  fact  that  a 
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vortex  at  a  downstream  station  of  the  apex  of  the  model  (ogive  cylinder)  will  react  to  the 
changed  flow  conditions  at  the  apex  at  a  time  At  later.  He  also  approximated  that  the 
disturbances  originating  at  the  apex  of  the  model  would  propagate  downstream  with  the 
free-stream  velocity  which  would  give  At  =  x/U^c-  Using  equations  6.11  and  6.12, 
the  effective  angle  of  attack  at  the  apex  of  the  model  {x  =  0.0)  can  be  obtained  as: 


* 

^effA  ^  ^  TI 

^oo 


(6.13) 


According  to  Ericsson’s  theory,  the  effective  angle  of  attack  aeff{x,t)  determining  the 
vortex  characteristics  at  station  x  downstream  of  the  apex,  is 


oteff{x,t')  =  aeffA{0,t-  At) 


(6.14) 


where  At  =  xfUoo 
attack  at  the  apex 
one  can  write 


By  using  the  Taylor  series  expansion  at  t  for  the  effective  angle  of 
and  assuming  small  non-dimensional  rotation  rates,  {aLlUoo)  <<  li 

=  (6.15) 


which  gives 


Oleff 


a  —  a 


Xcg  d”  X^ 


(6.16) 


Using  the  definition  given  by  the  equation  6.11,  we  can  obtain  Aae//  from  equation  6.16 


Aa.„  -  6,^^  (6.17) 

Equation  6.17  did  approximate  the  lags  in  asymmetric  vortex  shedding  on  the  pitching 
Ogive  cylinder  successfully. 

In  this  study,  Aae//  has  been  calculated  both  from  the  barebody  experimental  results  and 
equation  6.17  for  four  different  xjL  locations.  Since  the  pitch  rate  has  a  constant  value 
of  78  7s,  the  lag  equation  6.17  gave  constant  Aa*//  values  for  each  xjL  locations.  For 
calculating  Aae//  from  the  experimental  results,  certain  number  of  unsteady  separation 
locations  4>sep  and  corresponding  instantaneous  pitch  angles  a(t')  have  been  selected  from 
the  data.  Then  the  same  (l>sep  values  for  the  steady  data  and  the  corresponding  steady 
pitch  angles  have  been  picked.  For  a  specific  (f>sepi  fhe  difference  between  the  instanta¬ 
neous  pitch  angle  and  the  steady  pitch  angle  has  been  calculated  as  the  effective  angle  of 
attack  increment.  The  Aae//  results  obtained  from  equation  6.17  and  the  experimental 
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data  have  been  compared  in  figure  6.67.  As  can  be  seen  from  this  figure,  there  is  a  signif¬ 
icant  difference  between  experimental  Aoe//  and  the  one  calculated  from  equation  6.17. 
Especially  for  the  stations  xjL  >  0.501,  the  magnitudes  of  the  experimental  Aoe//  are 
twice  as  much  higher  than  that  of  obtained  by  using  the  model  equation  6.17.  Also, 
although  the  pitch  rate  is  constant,  experimental  Aoe//  does  change  with  the  instanta¬ 
neous  pitch  angle  thus  with  the  time.  This  comparison  shows  the  difference  between  the 
lags  in  the  unsteady  crossflow  separation  over  the  model  used  in  the  present  experiment 
and  the  flowfield  associated  with  the  asymmetric  vortex  shedding  on  the  ogive  cylinder. 
A  more  sophisticated  time-lag  modeling  technique  should  be  used  in  order  resolve  the 
complex  nature  of  the  time  history  effects  on  the  unsteady  cross  flow  separation  locations. 

6.3.2  First-Order  Differential  Time-Lag  Model 

Goman  and  Khrabrov  [35]  have  developed  a  first  order  time  lag  model  in  order  to  ap¬ 
proximate  the  time  history  of  a  dominant  flow  feature  such  as  the  separation  location 
in  general  unsteady  flows.  They  applied  this  model  to  pitching  two-dimensional  airfoils, 
delta  wings  and  the  unsteady  aerodynamics  of  a  complex  fighter  aircraft  configuration 
as  well.  They  used  the  separation  location  as  an  internal  state  variable,  and  defined 
the  forces  and  moments  as  functions  of  this  state  variable.  For  a  given  maneuver,  by 
obtaining  the  time  history  of  this  state  variable,  they  were  able  to  calculate  the  unsteady 
force  and  moments.  This  model  mainly  introduces  a  first-order  delay  differential  equa¬ 
tion  for  an  internal  state  variable  x  which  accounts  for  unsteady  effects  associated  with 
separated  and  vortex  flow.  The  variable  x  may  represent  the  location  of  trailing  edge 
flow  separation  on  a  pitching  airfoil  or  that  of  vortex  breakdown  on  a  delta  wing.  The 
original  form  of  the  differential  equation  governing  x  is  given  by; 

+  X  =  XQ{a  -  T2^)  (6.18) 

at  at 

where  ti  and  T2  are  time  constants  and  xq  describes  the  steady  dependency  of  x  on  a. 
Goman  and  Khrabrov  [35]  defined  as  a  quasi-steady  time  delay  associated  with  changes 
in  circulation,  boundary  layer  convection  lags,  and/or  boundary  layer  improvement  ef¬ 
fects,  while  Ti  as  a  relaxation  time  constant  associated  with  global  transient  aerodynamic 
effects,  or  dynamic  properties  of  the  separated  flow  adjustment. 

Wetzel  and  Simpson  [3]  implemented  an  extended  version  of  this  model  for  approximating 
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the  time  varying  nature  of  the  unsteady  separation  locations  over  the  maneuvering  pro¬ 
late  spheroid.  The  approximation  successfully  fit  the  experimental  unsteady  separation 
locations  and  non-dimensional  time  lag  values  for  different  x/L  stations  of  the  prolate 
spheroid.  The  same  extended  version  of  the  first-order  differential  lag  model  has  been 
used  in  this  study  in  order  to  approximate  the  unsteady  primary  separation  locations  of 
the  barebody  and  the  non-sail  region  of  the  sail-on-side  case: 

In  equation  6.19,  </>uns  represents  the  approximation  to  the  unsteady  separation  location 
and  t'  stands  for  the  first-order  non-dimensional  time  lag.  The  ^>0  is  the  quasi-steady 
separation  location  distribution  which  can  be  obtained  from  the  steady  separation  data 
at  each  a  =  o;(t').  Note  that  both  the  quasi-steady  separation  location  and  the  ap¬ 
proximation  to  the  unsteady  separation  location  do  also  vary  with  x/L.  This  is 
one  of  the  differences  between  the  extended  version  and  the  original  time  lag  model  of 
Goman  k.  Khrabrov  [35].  In  the  original  version,  a  single  point  of  separation  was  consid¬ 
ered.  However,  for  the  present  study  and  the  prolate  spheroid  case,  since  the  crossflow 
separation  occurs  along  a  line  rather  than  a  point,  (puns  and  (f>Q  were  also  functions  of 
x/L.  It  should  also  be  noted  that  the  time  lag  r'  taken  as  an  unknown  in  equation  6.19 
is  let  to  vary  in  x/L  and  identified  by  fitting  the  model  equation  with  the  experimental 
data. 


The  first-order  differential  equation  6.19  mainly  correlates  the  unsteady  separation  loca¬ 
tions  to  the  quasi-steady  data  by  the  time  lag  r'.  Note  that  if  we  solve  equation  6.19 
at  a  constant  x/L  location,  the  only  dependent  variable  would  be  the  non-dimensional 
time  t'.  Then  the  solution  to  this  initial  value  problem  can  be  obtained  by  first  taking 
the  Laplace  transform  of  each  side: 


^  */*uns(0)]  d”  0iiTts(s)  ^o(5) 


(6.20) 


After  rearranging  the  above  equation,  we  obtain: 

J.  („\  </’uns(0)  1  f  (po{s)  01\ 

Now  if  we  take  the  inverse  Laplace  transform  of  each  side  of  equation  6.21,  we  can  get 
the  final  expression  for  <j)uns{t')  as: 

k.n.(0)  +  i  r 

Jo 


4^uns{j'  ) 


(6.22) 


6.3  Time-Lag  Models 


135 


In  the  solution  procedure,  the  quasi-steady  separation  distribution  (po  was  obtained  by 
fitting  a  cubic  spline  to  the  steady  data.  The  independent  variable  was  t'  in  the  fitting 
procedure.  As  the  initial  value  for  the  first  steady  separation  location  in  the  steady 
data  was  used.  Since  equation  6.22  represents  the  solution  to  an  initial  value  problem,  the 
determination  of  the  initial  value  for  (puns  and  thus  the  first  steady  separation  location  is 
important.  However,  in  terms  of  getting  the  right  time  lag  value,  the  choice  of  the  initial 
value  does  not  have  a  crucial  role.  As  long  as  the  pitch  rate  a  is  constant,  the  same  time 
lag  will  be  obtained  even  if  the  initial  value  is  slightly  different  than  the  true  value  and 
the  uncertainty  in  the  initial  value  will  only  effect  the  initial  transient  part.  The  integral 
on  the  right  hand  side  of  the  equation  6.22  was  evaluated  numerically  for  each  t'  value 
by  using  the  trapezoidal  approximation.  For  each  x/L  station,  the  model  equation  6.22 
was  solved  with  different  values  of  the  r'  in  an  iterative  manner.  The  root  mean  square 
error  between  the  approximated  and  experimental  unsteady  data  was  calculated  at  each 
iteration  and  the  t'  that  gave  the  smallest  error  had  been  chosen  for  that  specific  xjL 
location. 

The  solid  lines  in  figures  6.40  through  6.59  show  the  results  of  the  time  lag  model  ap¬ 
proximation  to  the  unsteady  data  at  ten  x/L  locations  for  the  barebody  and  the  non-sail 
region  of  the  sail-on-side  case.  Note  that  all  the  stations  were  not  used  in  calculating  the 
time  lags.  At  the  stations  x/L  <  0.266,  no  unsteady  separation  location  at  any  instant  of 
time  was  observed.  At  the  last  station  x/L  =  0.863,  the  determination  of  the  unsteady 
separation  location  was  not  accurate  enough  because  of  the  flat  nature  of  the  Cf  profile 
near  the  minimum.  For  the  rest  of  the  stations,  the  first-order  lag  model  approxima¬ 
tion  fits  the  measured  unsteady  separation  locations  reasonably  well.  In  these  graphs, 
for  some  of  the  x/L  locations,  the  last  value  of  the  unsteady  and  the  steady  separation 
locations  are  slightly  different  although  they  should  be  the  same  ideally.  However,  the 
differences  are  within  the  uncertainty  of  separation  locations  that  is  ±2  °.  The  model 
fit  results  show  that  the  first-order  differential  lag  model  of  Goman  and  Khrabrov  does 
capture  the  time-varying  nature  of  the  unsteady  crossflow  separation  locations  reason¬ 
ably  and  may  be  used  as  part  of  the  unsteady  aerodynamic  models  used  to  describe  the 
physics  of  such  flows. 

Since  the  unsteady  separation  locations  do  not  follow  the  quasi-steady  separation  data 
with  a  time  lag,  a  first-order  lag  model  approximation  to  the  unsteady  data  of  the  sail 
side  cannot  be  obtained.  The  physical  model  given  by  the  differential  equation  6.19  does 
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not  describe  the  unsteady  separation  phenomena  of  the  sail  side.  In  order  to  approximate 
the  sail  side  unsteady  separation  pattern,  a  more  complicated  physical  model  should  be 
developed. 


Time  lags 

The  time  lag  values  r'  obtained  from  the  first-order  differential  lag  model  equation  fits  for 
different  xjL  locations  are  shown  in  figure  6.68.  In  this  figure,  circles  represent  the  time 
lags  of  the  barebody  case,  and  the  triangles  stand  for  the  results  of  the  non-sail  region  of 
the  sail-on-side  case.  For  the  barebody  case,  near  the  model  rotation  point  XcgjL  —  0.24, 
the  time  lag  value  is  close  to  zero  indicating  that  the  quasi-steady  data  follows  the 
unsteady  data  closely.  Then  an  increase  until  x/L  =  0.434  can  be  seen.  Since  there  are 
not  enough  points  in  this  region,  the  nature  of  this  increase  (linear  or  non-linear)  cannot 
be  determined  accurately.  After  this  point,  the  time  lag  stays  approximately  constant 
between  x/L  =  0.434  and  x/L  =  0.774,  taking  an  average  value  of  1.40.  At  x/L  =  0.819, 
this  value  drops  to  1.24.  The  non-dimensional  time  lag  t'  vs.  x/L  distribution  of  the 
non-sail  region  shows  a  similar  trend  as  the  one  obtained  for  the  barebody  case.  However, 
at  the  first  two  stations,  the  time  lag  values  are  considerable  higher  than  the  barebody 
values.  For  the  other  stations,  time  lags  are  approximately  constant  with  an  average 
value  of  1.50. 

Both  the  barebody  and  the  non-sail  region  t'  vs  x/L  distributions  are  different  fi:om  the 
one  obtained  for  the  prolate  spheroid  by  Wetzel  and  Simpson  [3].  In  that  case,  t'  changed 
approximately  in  a  linear  manner  increasing  in  downstream  direction  over  the  prolate 
spheroid  model.  By  making  a  linear  fit  to  the  time  lag  data  and  using  the  Ericsson’s 
Aaeff  approach  (see  equation  6.17)  with  an  unknown  convection  velocity  14  for  the 
propagation  speed  of  the  disturbances,  Wetzel  and  Simpson  [3]  were  able  to  determine 
that  the  disturbances  propagate  with  a  convection  velocity  of  Uc  —  0.3  x  Uoo-  It  should 
be  noted  that,  for  the  prolate  spheroid  pitch-up  maneuvers,  the  model  center  of  rotation 
point  was  at  x/L  =  0.5,  whereas  in  this  study  the  rotation  point  is  at  x/L  =  0.24. 
Also  the  Suboff  model  has  a  different  geometry  than  the  prolate  spheroid,  having  a  large 
proportion  of  a  constant  diameter  region.  Therefore,  one  may  expect  an  influence  of  the 
rotation  point  and  the  model  geometry  on  the  time  lags.  However  further  study  has  to 
be  done  in  order  investigate  the  true  complex  nature  of  the  time  lags  associated  with  the 
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unsteady  flow  field  over  this  undersea  vehicle  geometry. 
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Figure  6.1:  Comparison  of  barebody  Cf  vs.  (p  distribution  for  steady  and  unsteady  data  at 
all  x/L  locations,  a  =  3.1  t'  =  3.417,  a{t')  =  3.1  °. 
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Figure  6.2;  Comparison  of  barebody  C/  vs.  0  distribution  for  steady  and  unsteady  data  at 
all  x/L  locations,  o;  =  5.1  t'  =  3.818,  a{t')  —  5.1  °. 
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Figure  6.3:  Comparison  of  barebody  C/  vs.  distribution  for  steady  and  unsteady  data  at 
all  x/L  locations,  a  =  7.2  t'  =  4.238,  a{t')  =  7.2  ° . 
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Figure  6.4:  Comparison  of  barebody  C/  vs.  (p  distribution  for  steady  and  unsteady  data  at 
all  xjL  locations,  a  =  9.3° ,  t'  =  4.696,  a{t')  =  9.3 °. 
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Figure  6.5:  Comparison  of  barebody  C/  vs.  (j)  distribution  for  steady  and  unsteady  data  at 
all  xlL  locations,  a  =  11.3  t'  =  5.307,  a{t')  =  11.3  °. 
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Figure  6.7;  Comparison  of  barebody  Cf  vs.  4>  distribution  for  steady  and  unsteady  data  at 
all  xjL  locations,  a  =  15.3°,  t'  =  6.128,  a{t')  =  15.3°. 
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Figure  6.8:  Comparison  of  barebody  Cf  vs.  <j)  distribution  for  steady  and  unsteady  data  at 
all  x/L  locations,  a  =  17.4°,  t'  =  6.472,  a{t')  =  17.4°. 
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Figure  6.9:  Comparison  of  barebody  Cf  vs.  <j)  distribution  for  steady  and  unsteady  data  at 
all  x/L  locations,  a  —  19.4°,  t'  =  6.815,  ait')  =  19.4°. 
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Figure  6.10:  Comparison  of  barebody  Cf  vs.  0  distribution  for  steady  and  unsteady  data  at 
all  x/L  locations,  a  =  21.4°,  t'  =  7.235,  a{t')  =  21.4°. 
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Figure  6.13:  Comparison  of  barebody  C/  vs.  (/>  distribution  for  steady  and  unsteady  data  at 
all  xjL  locations,  a  =  27.6°,  t'  =  10.290,  oi{t')  —  27.6°. 
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Figure  6.14:  Comparison  of  sail-on-side  (region  without  the  sail)  Cj  vs.  ((>  distribution  for 
steady  and  unsteady  data  at  all  x/L  locations,  a  =  3.1  t'  =  3.293,  a{t')  =  3.1  °. 
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Figure  6.15:  Comparison  of  sail-on-side  (region  without  the  sail)  C/  vs.  <f>  distribution  for 
steady  and  unsteady  data  at  all  xjL  locations,  a  —  5.1  °,  t'  =  3.742,  a{t')  —  5.1  °. 
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Figure  6.16;  Comparison  of  sail-on-side  (region  without  the  sail)  C/  vs.  (j)  distribution  for 
steady  and  unsteady  data  at  all  xjL  locations.  a  =  7.2°,  t  =  4.261,  a{t )  —  7.2  . 
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Figure  6.17:  Comparison  of  sail-on-side  (region  without  the  sail)  C/  vs.  (j)  distribution  for 
steady  and  unsteady  data  at  all  x/L  locations,  a  =  9.3  t'  =  4.7123,  a{t')  =  9.3  °. 
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Figure  6.18:  Comparison  of  sail-on-side  (region  without  the  sail)  C/  vs.  (f)  distribution  for 
steady  and  unsteady  data  at  all  xjL  locations,  a  =  11.3  t'  =  5.197,  ot{t')  =  11.3°. 
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Figure  6.22:  Comparison  of  sail-on-side  (region  without  the  sail)  C/  vs.  4>  distribution  for 
steady  and  unsteady  data  at  all  x/L  locations,  a  =  19.4°,  t'  =  6.827,  =  19.4°. 
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Figure  6.25:  Comparison  of  sail-on-side  (region  without  the  sail)  Cf  vs.  (j)  distribution  for 
steady  and  unsteady  data  at  all  xjL  locations,  a  =  25.5  t'  =  8.266,  a{t')  =  25.5  °. 
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Figure  6.27:  Comparison  of  sail-on-side  (region  with  the  sail)  Cj  vs.  (j)  distribution  for  steady 
and  unsteady  data  at  all  xfL  locations,  a  =  3.1  t'  =  3.293,  a{t')  =  3.1  °. 
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Figure  6.29;  Comparison  of  sail-on-side  (region  with  the  sail)  Cf  vs.  4>  distribution  for  steady 
and  unsteady  data  at  all  x/L  locations,  a  =  7.2°,  t'  =  4.261,  a{t )  =  7.2  . 


Chapter  6.  Figures 


Chapter  6.  Figures 


r=5.197  arr)=7/.3  1 


steady 

unsteady 


3.(X)E-03 


iiooooooooooooooooooooooooooo89»«88§°°°°°°gS2^®®*** 

X/L=0.863  O.OOE+00  . ^ 

&oo o ooo 00000000000*888 8 8 88§»«9Qgooooooooo 8® 


X/L=0.819  O.OOE+00 


X/L=0.774  O.OOE+00 


X/L=0.706  O.OOE+00 


X/L=0.638  O.OOE+00 


X/L=0,570  O.OOE+00 


X  /L- 0.501  0. 00E-¥  00 


ooooooooooo 


ioooooSB9088889»#9qooooooooo28 


,.se88SSSSc 


TOOOOOOOOOOOOOOOOO 


. . 

I  o 888 o o oo8 80 •»8gooooog#§ 8 ooooooooooo 


POOOOOOOOOO 


Dooo88»»»«8oo°°°o88e»B80000g9®" 


f OOOOOOOOOOO 


o 


O  O  O  O 


O  o  O  O 


o  o  o  O 


O  O  O  o  O 


o  O  O 


Cl  g888888e890#«»*aaQgoo  g,ae8ooooo 

|.  J>oSS88888880«.»iQgo2S8* 


X/L=0.434  O.OOE^OO 


X/L^O.345  O.OOE+OO 


X/L^0.306  O.OOE^OO 


X/L=0.266  O.OOE-¥00 


X/L^0.170  Q.OOE^-00 


i888888e8se00**»«*e98ee 


.••  ••  ••• 

..•••••e®***®****®* 


—•••a  —qoC®® ®8* 

■»»e98B8999988»»®®^°*®*2  °  • 

*•  • 


o  O  o  o 


I  O  5 


•••*  ***8oO°° 

® 8800 000000°° 


*08 8 6 008 88888 8 8888 8888 8 8 000* 0000 00 0000 0000 00 0000 


8  8  8  8  0 


.*0e8888888SSS8SS58888888888888888888888888888888S  8  8  0  0  •  • 


X/L=0.I31  O.OOE+00 


X/L=0.110  O.OOE+00 


o80„,**08888888SSSSSSSSSSSSSSSSSSSSSSSSSSSSSoSSS;5SSSooS  88800. 

(  I  .  ■  I  ■  ■  -  .—t — ^ — I — ^ ^ i ^ ^ — > — ■ — ■ — ^ ^ — • — - - — ^ ^ ^ — 

180  200  220  240  260  o  2S0  300  320  340  360 

<t> 


Figure  6.31:  Comparison  of  sail-on-side  (region  with  the  sail)  C/  vs.  (j)  distribution  for  steady 
and  unsteady  data  at  all  x/L  locations,  oc  =  11.3°,  t'  =  5.197,  ot{t)  =  11.3  . 
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Figure  6.32:  Comparison  of  sail-on-side  (region  with  the  sail)  Cf  vs.  (j)  distribution  for  steady 
and  unsteady  data  at  ail  xjL  locations,  a  =  13.2  t'  =  5.663,  a{t')  =  13.2  °. 
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Figure  6.33:  Comparison  of  sail-on-side  (region  with  the  sail)  Cf  vs.  (j)  distribution  for  steady 
and  unsteady  data  at  all  rc/L  locations,  a  =  15.3°,  t'  =  6.070,  a{t')  =  15.3°. 
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Figure  6.36:  Comparison  of  sail-on-side  (region  with  the  sail)  Cf  vs.  4>  distribution  for  steady 
and  unsteady  data  at  all  x/L  locations,  a  =  21.4°,  t'  =  7.215,  a{t')  =  21.4°. 
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Figure  6.37:  Comparison  of  sail-on-side  (region  with  the  sail)  C/  vs.  0  distribution  for  steady 
and  unsteady  data  at  all  x/L  locations,  a  =  23.4  t'  =  7.626,  OL{i!)  =  23.4  °. 
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Figure  6.40;  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
xfL  =  0.266  for  the  barebody  case. 
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Figure  6.41:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.306  for  the  barebody  case. 
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Figure  6.42;  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.345  for  the  barebody  case. 
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Figure  6.43:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.434  for  the  barebody  case. 
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Figure  6.44:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.501  for  the  barebody  case. 
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Figure  6.45;  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.570  for  the  barebody  case. 
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Figure  6.46:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.638  for  the  barebody  case. 


Figure  6.47:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.706  for  the  barebody  case. 
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Figure  6.48:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.774  for  the  barebody  case. 


Figure  6.49:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.819  for  the  barebody  case. 
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Figure  6.50:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.266  for  the  non-sail  region  of  the  sail-on-side  case. 
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Figure  6.51:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.306  for  the  non-sail  region  of  the  sail-on-side  case. 


Figure  6.52;  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.345  for  the  non-sail  region  of  the  sail-on-side  case. 


Figure  6.53:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.434  for  the  non-sail  region  of  the  sail-on-side  case. 
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Figure  6.54;  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
xjL  =  0.501  for  the  non-sail  region  of  the  sail-on-side  case. 


Figure  6.55:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
z/L  =  0.570  for  the  non-sail  region  of  the  sail-on-side  case. 
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Figure  6.56:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
xlL  =  0.638  for  the  non-sail  region  of  the  sail-on-side  case. 


Figure  6.57:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.706  for  the  non-sail  region  of  the  sail-on-side  case. 


Figure  6.58:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.774  for  the  non-sail  region  of  the  sail-on-side  case. 


Figure  6.59:  First-order  differential  lag  approximation  to  the  unsteady  separation  data  at 
x/L  =  0.819  for  the  non-sail  region  of  the  sail-on-side  case. 
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Figure  6.60:  Steady  and  unsteady  separation  locations  (second  minima  in  Cj  measured  from 
(j)  =  180  °)  vs  t'  at  x/L  =  0.434  for  the  sail  region  of  the  sail-on-side  case. 
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Figure  6.61:  Steady  and  unsteady  separation  locations  (second  minima  in  C/  measured  from 
0  =  180  °)  vs  t'  at  x/L  =  0.501  for  the  sail  region  of  the  sail-on-side  case. 
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Figure  6.62:  Steady  and  unsteady  separation  locations  (second  minima  in  Cf  measured  from 
4>  =  180  °)  vs  t'  at  xlL  =  0.570  for  the  sail  region  of  the  sail-on-side  case. 


Figure  6.63:  Steady  and  unsteady  separation  locations  (second  minima  in  Cf  measured  from 
(f)  =  180  °)  vs  t'  at  x/L  =  0.638  for  the  sail  region  of  the  sail-on-side  case. 
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Figure  6.64;  Steady  and  unsteady  separation  locations  (second  minima  in  Cj  measured  from 
(f,  =  180  °)  vs  t'  dXxlL  =  0.706  for  the  sail  region  of  the  sail-on-side  case. 


Figure  6.65;  Steady  and  unsteady  separation  locations  (second  minima  in  C/  measured  from 
0  =  180  °)  vs  t'  at  x/L  =  0.774  for  the  sail  region  of  the  sail-on-side  case. 
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Figure  6.66:  Steady  and  unsteady  separation  locations  (second  minima  in  (7/  measured  from 
(f)  =  180  °)  vs  t'  at  x/L  =  0.819  for  the  sail  region  of  the  sail-on-side  case. 


Figure  6.67:  Incremental  effective  angle  of  attack  Aaeff  for  the  barebody  pitch-up  maneuver 
at  four  x/L  stations  as  a  function  of  instantaneous  angle  of  attack  Open  symbols  show 
the  Atte//  given  by  the  equation  6.17. (Both  angles  are  in  degrees) 
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Steady  and  the  unsteady  turbulent  surface  flow  on  a  generic  submarine  model  was  stud¬ 
ied.  Hot-film  sensors  with  constant  temperature  anemometers  were  used  to  measure 
the  steady  and  the  unsteady  wall  shear  magnitudes  over  the  model  surface.  Three- 
dimensional  separation  locations  were  determined  from  the  minima  of  the  skin-friction 
magnitudes.  Steady  skin-friction  magnitudes  were  obtained  at  fourteen  steady  angles 
of  attack.  The  dynamic  plunge-pitch-roll  model  mount  (DyPPiR)  was  used  to  simulate 
the  pitchup  maneuvers.  The  pitchup  maneuver  was  a  linear  ramp  from  1  °  to  27  °  in 
0.33  seconds.  Mean  wall  static  pressures  were  measured  at  10  °  and  20  °  angles  of  at¬ 
tack.  Qualitative  examination  of  the  steady  surface  flow  topology  was  done  by  using  the 
oil-flow  visualization  pictures. 

Two  model  configurations  were  studied  in  both  the  steady  and  the  unsteady  experiments: 
The  barebody  case  and  the  sail-on-side  case.  The  barebody  case  is  the  axisymmetric 
configuration  of  the  model  and  can  also  be  thought  as  a  missile  body  or  to  a  certain 
extent  fuselage  of  an  aircraft.  Therefore,  results  obtained  for  the  barebody  case  may 
be  applied  to  many  geometries  with  similar  shapes.  The  sail-on-side  measurements  were 
performed  with  the  sail  fixed  at  the  circumferential  position  (j)  =  270°.  Pitchup  maneuvers 
for  the  sail-on-side  case  simulate  the  turning  maneuver  of  a  submarine.  All  the  tests  were 
conducted  at  Rei  =  5.5  x  10®  with  a  nominal  wind  tunnel  speed  of  42.7±1%  m/s.  For  the 
experiments,  slotted  walls  with  an  open-air-ratio  of  38%  are  used  to  reduce  the  blockage 
effects.  However,  some  steady  pressure  and  hot-film  data  at  certain  angles  of  attack  were 
also  acquired  with  solid  walls  for  the  purpose  of  comparison. 
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For  the  range  of  conditions  at  which  the  tests  are  conducted,  steady  results  over  the 
barebody  show  that  the  flow  on  the  leeward  side  of  the  model  can  be  characterized  by 
the  crossflow  separation.  The  first  clear  detection  of  the  primary  separation  locations 
are  at  the  stations  downstream  of  xjL  =  0.638  at  a  =  5.1  As  the  angle  of  attack 
increases,  the  primary  separation  line  moves  upstream  and  at  a  specific  xjL  location 
moves  towards  the  windward  side.  At  a  =  11.3  °,  the  onset  of  the  secondary  separation 
can  be  observed  on  the  stern  region.  The  secondary  separation  line  migrates  upstream 
and  towards  the  leeward  side  with  increasing  angle  of  attack.  The  flat  circumferential 
skin-friction  profile  on  the  stern  region  indicates  the  low  speed  separated  flow  which  is  also 
consisted  with  the  oil-flow  pattern  of  this  region.  This  flat  profile  makes  the  identification 
of  the  true  minimum  and  thus  the  separation  location  difficult  for  the  stern  region.  A 
weak  separation  and  reattachment  of  the  flow  on  the  nose  region  of  the  model  can  also 
be  observed  at  high  angles  of  attack. 

The  results  of  the  sail-on-side  case  are  evaluated  in  two  separate  regions,  the  region  with 
no  sail  (between  0  =  0  °  and  180  °)  and  the  region  with  the  sail  (between  ^  =  180  °  and 
360°).  In  the  first  region,  the  origin  and  the  variation  of  the  primary  and  the  secondary 
separation  lines  as  a  function  of  xjL  and  a  show  the  same  characteristics  as  defined  for 
the  barebody  case.  This  implies  that  the  main  flow  feature  on  the  non-sail  region  is  the 
cross  flow  separation.  The  primary  separation  locations  of  the  non-sail  region  start  to 
deviate  from  the  barebody  results  having  an  offset  in  the  leeward  direction  at  angles  of 
attack  starting  from  a  =  15.3°.  On  the  other  hand,  secondary  separation  locations  of  the 
non-sail  region  are  shifted  in  the  leeward  direction  compared  to  the  secondary  separation 
locations  of  the  barebody  case  at  all  angles  of  attack. 

Downstream  of  the  sail,  the  flow  structure  on  the  sail  side  of  the  model  is  much  different 
than  the  one  observed  for  the  non-sail  side.  The  flow  field  in  this  region  is  strongly  affected 
by  the  presence  of  the  sail.  Compared  to  the  separation  topology  of  the  barebody  and 
the  non-sail  region  of  the  sail-on-side  case,  the  separation  location  trend  as  a  function 
oi  x/L  and  a  shows  significant  differences.  Two  minima  in  Cf  vs.  (j)  distributions  on 
the  leeward  side  of  the  sail  region  can  be  observed  for  certain  angles  of  attack  and  xfL 
locations.  The  first  minima  measured  from  (/>  =  180  °  can  be  observed  only  at  a  certain 
range  of  angles  of  attack,  while  second  minima  is  detected  at  all  angles  of  attack  starting 
from  a  =  5.1  °.  The  results  on  the  sail  side  indicate  that  the  flow  field  does  differ  from 
the  crossflow  separation  structure  observed  for  the  barebody  and  non-sail  region  of  the 
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sail-on-side  case.  Therefore  the  categorization  of  the  separation  locations  as  the  primary 
or  the  secondary  is  not  clear  and  may  not  reflect  the  real  flow  structure  of  this  region. 
The  flow  in  the  vicinity  of  the  sail-body  junction  is  dominated  by  the  horseshoe  type 
separation. 

The  comparison  of  the  barebody  pressure  measurements  made  with  the  solid  walls  and 
the  slotted  walls  shows  the  blockage  effect  in  the  solid  wall  case,  especially  at  a  =  20  °. 
On  the  leeward  side  of  the  model,  the  circumferential  pressure  gradient  is  found  to  be 
approximately  zero  over  the  regions  of  flow  separation.  One  can  think  of  using  the  flat 
pressure  profile  in  regions  of  separation  to  locate  the  separation  locations,  however  it  is 
difficult  to  determine  the  exact  point  where  this  flat  pressure  distribution  begins. 

Unsteady  skin-friction  distributions  of  the  barebody  case  and  the  non-sail  region  of  the 
sail-on-side  case  show  similar  trends  compared  to  the  steady  distributions  of  each.  For 
these  regions  at  high  instantaneous  angles  of  attack,  the  dominant  flow  feature  on  the  lee¬ 
ward  side  of  the  model  can  again  be  characterized  by  the  crossflow  separation.  However, 
the  unsteady  crossflow  topology  is  different  than  the  corresponding  quasi-steady  one;  the 
primary  separation  formation  occurs  at  higher  angles  of  attack  compared  to  the  steady 
ones  at  which  the  onset  of  the  first  steady  primary  separation  is  observed.  Also,  at  an 
instantaneous  pitch  angle  for  a  specific  xjL  location,  the  primary  unsteady  separation 
starts  more  leeward  compared  with  the  steady  case.  No  clear  secondary  separation  is 
observed  for  the  unsteady  maneuvers.  For  the  barebody  and  the  non-sail  region  of  the 
sail-on-side  case,  the  difference  in  the  separation  topology  originates  from  the  fact  that 
the  unsteady  separation  location  lags  the  unsteady  separation. 

As  also  seen  in  the  steady  case,  the  unsteady  separation  structure  of  the  sail  side  is 
different  from  the  unsteady  crossflow  separation  topology  observed  on  the  barebody  and 
the  non-sail  region  of  the  sail-on-side  case.  From  the  results  obtained  for  this  region,  it 
can  be  concluded  that  unsteady  separation  patterns  do  not  follow  the  quasi-steady  data 
with  a  time  lag.  A  lag  definition  as  described  for  the  unsteady  crossflow  separation  may 
not  be  appropriate  for  the  sail  side  separation  topology  based  on  the  results  obtained  in 
this  study. 

Time  lag  models  were  used  to  approximate  the  unsteady  separation  data  of  the  barebody 
case  and  the  non-sail  region  of  the  sail-on-side  case.  Algebraic  time  lag  model  results  do 
not  match  with  the  barebody  experimental  results.  This  indicates  the  complex  nature 
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of  the  unsteady  flow  separation  over  the  Suboff  model.  The  first-order  differential  time 
lag  model  of  Goman  and  Khrabrov  [35]  approximates  the  unsteady  data  reasonably  well 
and  captures  the  time-varying  nature  of  the  unsteady  crossflow  separation  locations. 
Therefore,  this  model  may  be  used  as  part  of  the  unsteady  aerodynamic  models  used  to 
describe  the  physics  of  such  flows. 

Time  lags  that  are  obtained  by  fitting  the  model  approximation  with  the  experimental 
unsteady  data  have  a  unique  variation  along  x/L  that  does  not  match  with  the  one 
obtained  from  the  prolate  spheroid  study  of  Wetzel  and  Simpson  [3] .  In  this  study,  time 
lag  values  are  approximately  constant  along  most  part  of  the  constant  diameter  region 
whereas  in  the  prolate  study,  a  linear  trend  increasing  in  the  downstream  direction  is 
observed.  However,  the  model  center  of  rotation  and  the  model  geometry  are  different 
between  two  cases.  This  may  imply  the  influence  of  these  two  parameters  on  the  time 
lags.  Further  study  has  to  be  done  in  order  to  investigate  the  effect  of  these  parameters. 
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Appendix  A 


Uncertainty  AncJysis  for  the 
Skin-friction  Measurements 


A.l  Classification  of  the  Uncertainties 

The  uncertainty  in  the  skin-friction  measurements  should  be  examined  in  two  levels.  The 
first  part  corresponds  to  the  uncertainty  in  the  calibration  of  the  hot-film  sensors  and 
can  be  thought  as  a  bias  source  of  uncertainty.  This  uncertainty  directly  affects  the  skin- 
friction  magnitudes,  however  it  does  not  have  an  effect  on  the  flow  separation  locations. 
The  second  one  is  the  random  relative  uncertainty  between  the  measurements  made  by 
one  sensor  at  different  circumferential  locations  on  the  model  in  different  flow  conditions. 
In  other  words,  this  uncertainty  determines  how  accurately  one  sensor  can  measure  the 
same  wall  shear  stress  value  in  a  given  circumferential  location  at  successive  measure¬ 
ments.  As  also  described  in  [2],  this  is  the  uncertainty  that  determines  the  accuracy  of 
the  flow  separation  locations. 

A. 2  Uncertainty  Calculations 

The  bias  (calibration)  uncertainty  is  mainly  related  to  the  accuracy  of  the  wall  shear  stress 
values  T-ui  used  in  the  calibration  process.  As  described  in  chapter  4,  these  shear  stress 
values  were  calculated  by  using  the  boundary  layer  properties  that  were  obtained  from 
the  boundary  layer  velocity  profile  measurements.  The  previous  experiments  that  used 
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the  same  technique  approximates  an  uncertainty  of  ±5%  in  the  skin  friction  coefficient 
Cf  (Simpson  [36]).  Since  the  uncertainty  in  the  determination  of  the  free-stream  velocity 
5Uoo  was  less  than  1%  during  the  measurements,  the  uncertainty  of  the  wall  shear  stress 
5tu)  can  also  be  thought  as  ±5%.  Thus,  this  value  may  be  taken  as  the  bias  uncertainty 
in  the  calibration  process. 

In  order  to  calculate  the  relative  uncertainty  in  the  actual  measurements,  we  should  again 
consider  the  King’s  equation  written  for  the  constant-temperature  hot-film  anemometers 
(CTA): 

rl/^  =  A^  +  B  (A.l) 

where  d  is  the  difference  between  the  sensor  operating  temperature  (which  is  constant 
for  CTA)  and  the  free-stream  temperature: 

d  =  n-T^  (A.2) 


The  relative  uncertainty  of  measured  by  each  sensor  will  be  due  to  the  random  un¬ 
certainties  of  the  measured  sensor  voltage  5E  and  the  temperature  difference  69.  Note 
that  the  calibration  constants  A  and  B  are  not  considered  here,  since  they  don’t  have  a 
random  variation.  These  coefficients  have  the  same  value  between  each  calibration  run 
as  described  in  chapter  4.  The  overall  relative  uncertainty  in  can  be  determined  by 
using  the  Kline  and  McClintock  method  given  in  Doebelin  [37]: 


In  the  above  equation,  W  =  W{xi,X2,  ■ . .  ,Xn)  represents  the  measured  variable;  5W, 
the  uncertainty  of  the  measured  variable;  Xj,  each  variable  that  effect  the  measured 
quantity;  and  5xi,  the  uncertainty  in  each  of  these  variables.  The  partial  derivatives  are 
the  sensitivities  of  W  with  respect  to  each  Xi.  In  order  to  use  this  method,  each  6xi 
should  be  random  and  uncorrelated  and  have  the  same  confidence  level  (e.g.  20:1  odds). 
Then  the  overall  random  uncertainty  in  W  will  also  have  the  same  odds. 


By  using  equation  A. 3,  can  be  approximated  by: 


(A.4) 


Since  the  sensor  temperature  is  constant,  temperature  difference  uncertainty  will  be 
equal  to  the  uncertainty  of  the  free-stream  temperature  which  was  at  most  ±0.5  °C 
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between  each  calibration  run  in  the  experiments.  By  using  the  nominal  1.1  value  of  the 
overheat  ratio, 

^  =  ±0.013  (A.5) 

6 

can  be  obtained.  In  order  to  get  statistically  meaningful  uncertainty  results,  equation  A.4 
has  been  calculated  for  each  steady  and  unsteady  skin-friction  value  in  the  corresponding 
steady  and  the  unsteady  data  reduction  programs.  The  value  of  each  5E‘^  has  been 
obtained  by  using  the  voltage  standard  deviation  info  calculated  in  the  same  programs. 
At  the  last  step,  the  overall  random  uncertainty  of  Tw  was  obtained  by  considering: 


STyj  _  5tJ 

J.  ~  1/3 

Till 


(A.6) 


Note  that  the  uncertainty  in  the  skin- friction  5Cf  is  equal  to  the  wall  shear  stress  un¬ 
certainty  St-uj,  since  the  uncertainty  of  the  free-stream  velocity  is  negligible.  As  the  final 
results,  the  random  uncertainties  in  Cf  with  20:1  odds  were  obtained  as: 


for  the  steady  measurements  6Cf  =  6%  of  C/ 
for  the  unsteady  measurements  6Cf  =  8%  of  Cf 


The  higher  uncertainty  in  the  unsteady  measurements  come  from  the  relatively  higher 
unsteady  voltage  uncertainties.  Above  values  lead  to  an  uncertainty  of  ±2  °  in  the 
determination  of  the  separation  locations. 
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Robust  Locally  Weighted  Regression 
and  Smoothing  LOESS 


Robust  locally  weighted  regression  technique  (Loess)  is  a  smoothing  method  which  uses 
the  pre-determined  windows  of  the  original  data  in  the  regression  process.  Within  each 
individual  window,  the  data  points  near  to  the  location  where  the  smoothed  value  of 
the  dependent  variable  is  sought  are  given  higher  weights.  Outliers  in  the  data  set  are 
detected  during  the  fitting  process  in  an  iterative  manner  and  not  used  in  the  regression. 
These  features  of  Loess  technique  enable  to  smooth  the  data  without  changing  the  original 
pattern. 

B.l  Mathematical  Description 

Mathematical  details  of  the  Loess  method  are  described  in  Cleveland  [29]  and  [30].  The 
theory  behind  the  smoothing  process  can  be  summarized  as  follows:  Let  {x,y)  be  the 
point  where  we  seek  for  the  smoothed  value  g{x).  Here  x  corresponds  to  the  independent 
and  y  to  the  dependent  variable  in  our  data  set.  Two  parameters  are  needed  to  be 
chosen  to  fit  a  loess  curve.  The  first  parameter,  a,  is  a  smoothing  parameter;  it  can 
be  any  positive  number  but  typical  values  are  0.25  to  1.0.  As  a  increases,  the  curve 
becomes  smoother.  The  second  parameter.  A,  is  the  degree  of  the  certain  polynomials 
that  are  fitted  by  the  method;  A  can  be  1  or  2.  Suppose  first  that  a  <  1.0.  Let  q  be  an 
truncated  to  an  integer  where  n  is  the  total  number  of  the  data  points.  We  assume  that 
a  is  large  enough  so  that  q  is  at  least  1,  although  q  is  much  larger  than  1  in  most  of  the 
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applications. 

Let  Ai(x)  =  |.T,-  -  .t|  be  the  distance  from  x  to  Xi,  and  A(i)(.T)  be  these  distances  ordered 
from  the  smallest  to  the  largest.  T(u),  the  tricube  weight  function  is  defined  as: 

(1  -  \u\^f  for  |ri|  <  1.0 

0.0  otherwise 


Then  the  neighborhood  weight  given  to  the  observation  {xi,yi)  for  the  fit  at  x  is 


Wi(x)  =  T 


(  ^i{x)  \ 


(B.2) 


For  Xi  such  that  Aj(a:)  <  A(,)(x),  the  weights  are  positive  and  decrease  as  Ai(x)  increases. 
For  Ai(x)  >  A(g)(x),  the  weights  are  zero.  If  a  >  1,  the  Wi{x)  are  defined  in  the  same 
manner,  but  A(g)(x)  is  replaced  by  A(„)(x)o!. 

If  A  =  1,  a  line  is  fitted  to  the  local  data  using  weighted  least  squares  with  weight  Wi{x) 
at  {xi,yi)\  and  the  values  of  a  and  6  are  found  that  minimize: 


^  Wi{x)  {yi-a-  bxif  (B.3) 

1=1 


Let  a  and  b  the  minimizing  values,  then  the  fit  at  x  is 


g{x)  =  a +  bx 


(B.4) 


For  A  =  2,  a  quadratic  polynomial  is  fitted  to  the  local  data  using  weighted  least  squares; 
values  of  a,  6,  and,  c  are  found  that  minimize 

Wi{x)  {yi  -  a-bxi-  cxf)^  (B.5) 

i=l 

If  a,  6,  and  c  are  the  minimizing  values,  then  the  fit  at  x  is 

g{x)  =  a  +  6x  +  cx^  (B.6) 


B.2  Selection  of  the  Loess  Parameters 

In  the  previous  section,  two  important  parameters  a  and  A  that  must  be  chosen  in  the 
smoothing  process  are  defined.  As  described  in  Cleveland  [29],  in  any  specific  application 
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of  the  loess  method,  the  choice  of  these  two  parameters  must  be  based  on  a  combination  of 
judgment  and  of  trial  and  error.  As  o  increases,  the  loess  fit  becomes  smoother.  However 
after  a  certain  value  of  o,  a  lack  of  fit  may  be  observed  as  the  peaks  and  valleys  in  the 
data  are  missed  because  of  the  large  value  of  cx.  For  small  values  of  a,  the  underlying 
pattern  is  tracked,  but  the  local  wiggles  do  not  appear  to  be  supported  by  the  data.  The 
goal  in  choosing  a  is  to  produce  a  fit  as  smooth  as  possible  without  unduly  distorting  the 
underlying  pattern  in  the  data.  Cleveland  [29]  also  shows  that  a  locally  quadratic  fit 
(A  =  2)  may  be  superior  to  the  linear  fit  (A  =  1)  in  terms  of  picking  the  local  minima 
and  maxima  in  certain  data  sets. 

In  this  study,  the  loess  method  is  used  in  smoothing  some  Cf  vs  (p  distributions,  and  the 
optimum  combination  of  the  parameters  for  each  configuration  and  test  condition  have 
been  found  to  be  as  follows: 

Unsteady,  barebody  a  =  0.25  and  A  =  1 

Steady,  sail-on-side  a  =  0.10  and  A  =  1 

Unsteady,  sail-on-side  (sail  region)  a  =  0.20  and  A  =  1 

Unsteady,  sail-on-side  (non-sail  region)  a  =  0.25  and  A  =  1 

Figures  B.l  and  B.2  show  the  loess  fits  to  the  barebody  unsteady  Cj  vs.  (p  data  at 
a{t')  =  11.3  °  and  x/L  =  0.819,  and  at  a{t')  =  21.5  °  and  xjL  =  0.638  respectively.  As 
can  be  seen  from  these  figures,  by  using  the  loess  parameters  A  =  1  and  a  =  0.25,  a 
desirable  level  of  smoothing  is  achieved  without  changing  the  general  trend  and  the  local 
minima  of  the  data. 


